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A detailed analysis of the growth of a Bose-Einstein condensate is given, based on quantum kinetic theory, 
in which we take account of the evolution of the occupations of lower trap levels, and of the full Bose-Einstein 
formula for the occupations of higher trap levels, as well as the Bose stimulated direct transfer of atoms to the 
condensate level introduced by Gardiner et al. We find good agreement with experiment at higher temperatures, 
but at lower temperatures the experimentally observed growth rate is somewhat more rapid. We also confirm 
the picture of the "kinetic" region of evolution, introduced by Kagan et al, for the time up to the initiation 
of the condensate. The behavior after initiation essentially follows our original growth equation, but with a 
substantially increased rate coefficient. 

Our modelling of growth implicitly gives a model of the spatial shape of the condensate vapor system as the 
condensate grows, and thus provides an alternative to the present phenomenological fitting procedure, based 
on the sum of a zero-chemical potential vapor and a Thomas-Fermi shaped condensate. Our method may give 
substantially different results for condensate numbers and temperatures obtained from phenomentological fits, 
and indicates the need for more systematic investigation of the growth dynamics of the condensate from a 
supersaturated vapor. 



I. INTRODUCTION 

Although the race to produce a Bose-Einstein condensate 
was preceded by intense debate concerning the likely rate of 
its formation, the discovery that a Bose-Einstein condensate 
of alkali atoms could be produced relatively simply [Q-gjl, and 
that the growth time was of the order of one second, moved 
most theoretical activity into the investigation of the proper- 
ties of the condensates so produced. Since the production of 
the first Bose-Einstein condensate there have been few the- 
oretical investigations into condensate growth, and only one 
experiment [Q] has made any measurements of growth rates. 
Only the work of the present authors and co-workers, based 
on quantum kinetic theory, has made quantitative predictions 
on the growth rate of a Bose-Einstein condensate. This work 
started when we showed how to introduce the concept of stim- 
ulated condensate growth resulting from kinetic processes [||], 
leading to a very simple formula for the growth rate. The MIT 
experiment [Q] took the form of a verification of the validity 
of our theoretical prediction. At the same time in [^] we re- 
fined the basic concept of Bosonic stimulation to generate a 
less idealized theoretical picture, and to compare it with ex- 
periment. These initial papers were of necessity brief, and 
developed neither the full theoretical justification on the nu- 
merical modeling nor the full range of possible comparison 
with the available experimental data. In particular, no account 
was taken of the information available on the spatial distribu- 
tion of the atoms in the vapor-condensate system as the con- 
densate grows from the vapor. 

This paper will therefore give the detailed justifications and 
a full range of comparison with experimental data. Most par- 
ticularly, we want to present a theoretically justifiable method 
of describing the condensate vapor system as it grows. The 
absence of such a description has led to a phenomenologi- 
cal fitting of vapor profiles to a zero chemical potential |^,^|] 
Bose-Einstein distributions, which may be an imperfect model 



whose results could well be misleading. 

The theoretical description of condensate growth that we 
present is largely able to be viewed as a modification of the 
quantum Boltzmann equation, in which, however, explicit 
note is taken of the modification of the excitation spectrum by 
the existence of the condensate, including of course the fact 
that the lowest single particle excitation energy is the chemi- 
cal potential pc{no) of the condensate of no atoms. Equilib- 
rium arises as a result of the equality of the chemical poten- 
tials of uncondensed vapor and condensate, a picture which 
is rather similar to that normally adopted for chemical reac- 
tions. The quantum Boltzmann equation itself automatically 
provides the Bose stimulation, which makes transition rates 
into the condensate and other highly occupied levels achieve 
a speed which permits the production of the condensate in 
a finite time. Without Bose stimulation, the production of a 
condensate of about 1,000,000 sodium atoms would take 30 
hours, rather than the 100ms observed. 

At first glance it might appear that a description which ap- 
pears to be based on the quantum Boltzmann equation would 
have nothing to say about condensate coherence or the ori- 
gin of that coherence. This is emphatically not the case — 
the kinetics of the transfer of the between energy levels in 
a trap requires the existence of a wavefunction for each en- 
ergy level. The condensate level has its own wavefunction, 
and this obeys the Gross-Pitaevskii equation. The coherence 
arises because this level becomes macroscopically occupied. 
There is no precise moment when one can say that the con- 
densate initiates. This picture applies in a trap, in which the 
energy levels about which we have been speaking are rather 
well separated. The picture of a Bose-Einstein condensate 
as developed in the middle part of this century as a part of 
condensed matter theory is of a homogeneous and thus in- 
finitely extended system — a system for which the thermody- 
namic limit is achieved. Looked at from our viewpoint, this 
would be achieved by making the trap broader and ultimately 
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flat. There is a transition point where the trap becomes so flat 
that there is an occupation of the lowest quasiparticle levels 
which becomes comparable with the occupation of the con- 
densate itself. At this stage the traditional condensed matter 
picture becomes relevant, but this is not achieved in any traps 
presently in use. 



II. MODEL FOR GROWTH OF A CONDENSATE 

In this section, the formalism of Quantum Kinetic theory 
|]|], will be used to form a model of the growth of a trapped 
Bose-Einstein condensate. The Bose atoms are described by a 
second-quantized field, in the pseudopotential approximation; 
that is, we write 
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and the term Ht arises from a trapping potential as 

H T = J d 3 xl/ T (x)V' t (x)V(x). 

The pseudopotential method is used — its validity for this kind 
of system has been justified in QKV — where u = Airah 2 /m, 
and a is the s-wave scattering length arising from the inter- 
atomic potential. 

The situation being considered is that of a vapor cloud con- 
fined in a trap in which the lower energy levels are not sig- 
nificantly populated, whilst the higher energy levels contain 
thermalized equilibrium populations, characterized by a tem- 
perature T and chemical potential fj,, unstable against conden- 
sate formation. 

This situation is likely to arise, to a degree of approxima- 
tion, if a system, which is initially in equilibrium at a temper- 
ature slightly greater than the critical temperature, is cooled 
very suddenly to a temperature below the critical temperature, 
by removing the very high energy atoms in a rapid evapora- 
tive cooling 'cut' . The higher energy levels will very quickly 
come to their equilibrium distributions, since the difference 
between the distributions before and after the cut are quite 
small at the higher energies. The lower levels will however be 
far from equilibrium and evolve to form a condensate. This 
is possibly one of the easiest scenarios to model, it is also the 
situation investigated by the only detailed experimental study 
of the growth of a condensate in a gas of 23 Na atoms [H], 

A. The System 

The system is contained in a three-dimensional harmonic 
potential, characterized by the frequencies co x , lo v and u) z . It 
will be useful to define the geometrical mean frequency of the 
trap as lj = (wjW^) 1 / 3 . 



1. Effective potentials arising from mean-field effects 



In this system the energies and wavefunctions of the lower 
trap levels are quite strongly affected by the presence of the 
condensate, and the effect will of course change as the con- 
densate grows. In QKV it was shown that it is reasonable 
to account for this by introducing mean-field effects, which 
make the effective potentials depend on the occupations of 
the bands. The situation is illustrated in Fig.|l} As the conden- 
sate grows, it expels the vapor from the center of the trap, and 
this expulsion serves to reduce the mean-field of the vapor as 
experienced by the condensate. The growth will be assumed 
to be so slow that the condensate- and non-condensate-bands 
are always in thermal equilibrium — that is, they will have a 
well defined temperature shared by both of them, but will not 
have the same chemical potential. Growth therefore occurs as 
atoms are transferred from the vapor to the condensate, lead- 
ing eventually to a unique chemical potential for the whole 
system. 
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FIG. 1. Representation of the modification of the trapping poten- 
tial for a) the non-condensate band and b) the condensate band due 
to mean-field effects. 



2. Condensate- and non-condensate bands 



In the formalism of quantum kinetic theory, the system is 
divided into condensate and non-condensate bands. We will 
treat in this paper the situation in which the non-condensate 
band is assumed to be in thermal equilibrium with a temper- 
ature T and chemical potential (j,, and to contain the vast ma- 
jority of the atoms so that it is essentially undepleted by the 
process of condensate formation. 

The picture of growth we will use is that presented in QKV. 
In that paper, it was shown that a legitimate division into con- 
densate and non-condensate bands can be made in which one 
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distinguishes between particle-like excitations, to which it is 
possible to assign a definite number of atoms, and phonon-like 
excitations, which are collective modes, which normally in- 
volve a large average number of atoms, but are not eigenstates 
of the atom number. In practice, it has been shown [|l4|] that 
the energy above which all excitations are essentially particle- 
like is relatively small. 

For the purposes of our modeling, however, there are two 
criteria which must be considered in the definition of the con- 
densate band. 

i) The non-condensate band is considered to be time inde- 
pendent, therefore the condensate band must include all levels 
whose populations change significantly during the condensate 
growth process. For the non-condensate band, the thermal 
distribution is given in the bulk by \e( E ~^/ kT — This 
is only valid for E > fj,, and gives very large populations 
when E w /i. The transition rates in and out of levels in 
this vicinity also become very large, which contradicts the as- 
sumption that the distribution of the non-condensate band is 
time-independent. These lower states are therefore must be 
treated time dependently, and hence must be included in the 
condensate band. 

ii) The condensate band consists principally of levels whose 
energy eigenvalues are significantly affected by the presence 
of a condensate — but levels which are not affected may be 
included if this is desirable, which must be done if the first 
criterion is to be met. 

Consequently, in this paper we will choose the condensate 
band to consist of all levels with energy less than the value 
Er. We will also introduce an energy E& < Er, which is 
the energy above which we can consider the energy levels to 
be unaffected by the condensate, as illustrated in Fig.|[ 



3. Grouping of energy levels into bands 

The inclusion of all the condensate band energy levels in 
the model means that simulations of the system require, in 
principle, the calculation of all the eigenfunctions of the con- 
densate band, and detailed summations over these. In prac- 
tice the number of energy levels involved is of the order of 
tens of thousands, which makes an exact description imprac- 
tical. However, progress can be made by grouping together 
energy levels in the condensate band into small 'sub-bands', 
with only the ground state (condensate state) being described 
as a single level. Each sub-band is described by an energy 
e m and contains all the eigenstates found within the energy 
range [e m — Ae TO /2, e m + Ae m /2]. The value of Ae m is cho- 
sen partially by the requirement that the lowest of these sub- 
bands contains at least 3 levels. Smaller values of Ae m would 
lead sub-bands containing only fractions of individual levels, 
which is obviously unphysical. As the condensate grows, the 
mean-field effects from the high occupation of the condensate 
level will cause the energies of the levels in the sub-bands to 
increase. The values of e m and Ae m are therefore dependent 
on the condensate occupation, and the manner in which they 
are altered will be discussed later. 



B. Notation 

For clarity, we set out some of our notation. 

TV : number of atoms in the condensate band, (5) 

no : number of atoms in the condensate, (6) 

no j : equilibrium number of atoms in the condensate, (7) 

Hc{no) '■ chemical potential of the condensate, (8) 

fj, : chemical potential of the non-condensate band, (9) 

£ no (x) : wavefunction of an n atom condensate. (10) 

In the situations we will consider, the number of condensate 
atoms no will vary from zero to almost TV, but this will always 
be substantially less than the number of atoms in the whole 
system, composed of both condensate and non-condensate 
bands. Thus, when the condensate is fully grown, the approx- 
imation no ~ TV will be valid, and will often be used. 



C. Density of States for the System 



In the absence of any condensate, the density of states 
G(E) is taken to be that of a non-interacting gas in a harmonic 
well. That is 



(E - 

2h 3 LU x LO y U! z 



(11) 



where N(75) is the cumulative number of states with en- 
ergy less than E and Q = (oj x + w y + u) z )/3. The num- 
ber of states in the sub-band with average energy e m is thus 
g m = G(e m )Ae m . The energy scale is such that the value of 
Vt{y) = m(u> 2 x 2 + u 2 y 2 + cu 2 z 2 )/2 is zero at the origin. 

Once the condensate begins to form, the mean field effects 
need to be taken into account. The mean field repulsion due 
to the condensate changes the energies of the lower trapped 
states. The energy of the condensate level is equal to the 
chemical potential uc(no), which increases with no. In the 

Thomas-Fermi approximation the rise is proportional to n^ 5 . 
Obviously the condensate level must remain the lowest en- 
ergy state, and thus the energies of the other states below E& 
must also rise in some fashion. The exact nature of the energy 
change is difficult to calculate, but some reasonable approxi- 
mations can be proposed which should reproduce the signifi- 
cant behavior caused by the mean field effects. 

The most simply calculated estimate of the energy changes 
is to assume the energies of the sub-bands e m are evenly dis- 
tributed between the fixed upper limit of Ea and the lower 



limit of uc (no) 



an 
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This is illustrated schematically in 



Fig.[| Both the values of e m and the values of Ae m are now 
no dependent. 

The final value of Ae m (after the growth of the condensate) 
is set to be fioj. This condition also always fulfills the require- 
ments that there are at least about 3 discrete levels contained in 
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the lowest energy sub-band, and yet ensures that the sub-band 
has only a relatively small energy range. 

The density of states for the condensate band, in the pres- 
ence of the condensate, is thus taken to be approximately 
G no (E) = JV[E — nc( n o)] 2 , where AT is a normalization 
chosen so that the cumulative number of states at Ea is the 
same as for the non-interacting harmonic oscillator potential. 
This behavior is illustrated in Fig.||. As the discontinuity at 
En shows, this model is obviously quite simplistic, and more 
realistic models will be discussed later. 

It should be noted that in the inset of Fig.[3] the number of 
particles per energy interval f(E) is shown (the occupation of 
the condensate level is not shown) for equilibrium conditions. 
From this inset it can be seen that the vast majority of atoms 
do indeed reside at energies higher than Er, and so the as- 
sumption that the non-condensate band is undepleted should 
be valid. 



I 

Non-Condensate 
Band 



Condensate Band 



FIG. 2. Schematic representation of the simple model to change 
the energy levels due to the mean field effects of the condensate. Lev- 
els are evenly distributed between Ea and (ic(no) which increases 
with no. 



1. Modified Thomas-Fermi Chemical Potential 

In the Thomas Fermi approximation the chemical potential 
of the condensate is given by 




Mc(no) 



1 5 auj x uj y lo z m 1 / 2 Ti 
2^2 



no 



2/5 
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which vanishes as no — > 0. However, the ground state of a 



non-interacting gas in a harmonic well is %(u 



-OA 



and so the real chemical potential should approach this value 
as no — > 0. In order to interpolate the Thomas-Fermi chem- 
ical potential to satisfy this requirement, the following form 
for the chemical potential will be used 



(13) 



where a = (lbaLUxLUyUJzm 1 / 2 ^ 2 /A^/2) 2 / 5 , and v is a con- 
stant such that av 2 l' J = %{lo x + u y + uj z )/2. 



2. Estimate of Ea 

An estimate for the value of Ea, above which the exci- 
tation spectrum is well described by that for a non-interacting 



gas, can be obtained using the number-conserving Bogoliubov 
spectrum [|l5]Jll]]. The quantity of interest now is the ratio of 
the corrections to the energy level arising from the presence 
of the condensate, to the energy level determined by the non- 
interacting gas model. For the case of the trap used in the 23 Na 
growth experiments at MIT [Q], numerical calculations show 
that the ratio is less than 10% for energies > 2/ic (no), and the 
corrections are of the order of only 5% for E ~ 2.5uc(no). 
Thus a reasonable estimate of Ea to be used in the simula- 
tions is 



Ea = (no,/) 



(14) 



where no./ is the equilibrium occupation of the condensate 
level, and this is the value of Ea which will be used in this 
paper. 
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FIG. 3. The cumulative number of states for a gas in a har- 
monic potential well. The dotted line shows the situation for 
a non-interacting gas (no mean field effects). The dashed line 
represents the cumulative number of states due to the simple 
model proposed in the text to incorporate mean field effects. The 
solid line shows the corresponding occupation per energy interval 
f(E) = G(E)[exp((E - (j)/kT) - also shown on a larger 

scale in the inset. 



3. Comparison with More Accurate Density of States 

As a result of the predominantly single particle nature of the 
excitation spectrum, the cumulative number of states N(i?) is 
expected to be quite well described at high energies by the 
semiclassical approximation 



N(E) = 



1 



(2tt^) : 



J de J dr J dp S(e - E sp (p, r)) 



(15) 



where E sp {p,r) = p 2 /2m + V T (r) + (8irh 2 a/m)\£ o (r)\ - 
Hc(n-o) is the semiclassical energy of a single particle in the 
potential created by the combination of the trapping poten- 
tial and the mean field repulsion of the condensate. Equa- 
tion ( |l5| ) thus represents the summation over all phase space 
cells which contain a single particle excitation state of energy 
less than E. In flip]], by carrying out the momentum and space 
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integrals, an expression for equation ( |l5| ) was found, for the 
case of an isotropic harmonic trap of frequency u>, in the form 

N(E) 



no 



E 2 d(l-x)(x + r]/E) 



1 , . [x 2 + e 2 / V 2 } 1 / 2 -x 
ax er]\ 



where 



kT r , ' 



( n a 



VC(3) 



■ E 2 jrj 1 



1/3 



(16) 



(17) 



The chemical potential (ic(n>o) is given by the Thomas-Fermi 
approximate form ( |l2| ) and the energy is given in the dimen- 
sionless units E = E/kT c . This semiclassical form for N(E) 
was found in [ [T6] , |l4] ] to be practically indistinguishable from 
that found by numerical solutions of the Bogoliubov spectrum 
over the entire range of energies. 
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Energy in units of p c (n Q ) A 
FIG. 4. The cumulative number of states N(_B) below an energy 
E. The solid line shows the results obtained using the semiclassi- 
cal approximation for an isotropic trap equation (jig), the dotted line 
shows the case for the non-interacting harmonic well. The dashed 
line represents the form used in this paper with Ea = 2pc(no)- 
The results were obtained for a condensate of 5 million atoms at a 
temperature of 900nK. 

In Fig.Hthe semiclassical form of N(E) obtained from 
equation ( jl6| ) is compared to that of the non-interacting har- 
monic oscillator (given by equation (pi])), and the density of 
states obtained using our model with Ea equal to 2(xc(no). 
The figure does not show very good agreement of our model 
with the semiclassical results at moderate energies, although 
at low energies the agreement is good. At high enough en- 
ergies (not shown) the non-interacting potential results be- 
come practically indistinguishable from those obtained from 
the semiclassical method, the energy at which this occurs is 
about 5fic( n o) for the results in Fig||. It should be em- 
phasized that this semiclassical result applies to an isotropic 



trap, and the consideration of the anisotropy of realistic traps, 
which has not yet been accounted for, may have a significant 
effect on the spectrum. 



4. Approximations 

Several approximations that have been made in the deriva- 
tion of this model rely on the condensate band being small 
relative to the non-condensate band, and if Er and Ea are 
too large then these approximations will not be valid. The 
approximations concerned are: 

i) That the non-condensate band is so large that it is essen- 
tially undepleted by the process of condensate growth. 

ii) That the scattering processes between condensate band 
atoms may be taken as being negligible compared to the inter- 
band scattering processes, as will be assumed later. 

It should be noted, and will be shown later, that the major 
effect on the overall growth due to mean field effects is caused 
by the changes in energies of the lowest energy levels, and for 
these levels the model proposed here is in quite good agree- 
ment with the semiclassical results. Because of this, as well 
as for the above reasons, the value of E& used will generally 
remain equal to 2fic(n j). The value of Er is chosen some- 
what larger. This provides a check that the solutions we find 
do match smoothly onto the distribution above En, which is 
assumed not to change. 
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FIG. 5. Illustration of the two types of dynamical processes, 
growth and scattering. Only the two non-condensate particle 
growth process is shown (process (a) in the text), and only 
non-condensate/condensate band scattering is shown (process (c)). 



D. Dynamical Processes 

The dynamics which will be considered in order to describe 
the evolution of the condensate band arise from the following 
processes: 
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a) Two particles in the non-condensate band collide, one of 
the particles leaves with an increased energy, and the remain- 
ing particle enters the condensate band, having now an en- 
ergy less than Ef> . Of course the reverse process must also be 
considered — a non-condensate band particle colliding with a 
condensate band particle and exciting it out of the condensate 
band. 

b) A non-condensate band particle collides with a condensate 
band particle and exchanges energy such that both particles 
end up in the condensate band, and the reverse process. 

c) A non-condensate band particle collides with a condensate 
band particle, transferring some energy, but both particles re- 
main in their respective bands. 

d) Two particles in the condensate band collide, transferring 
energy, with the result that both particles remain in the con- 
densate band, but having different energies than before the 
collision. 

Processes (c) and (d) will be termed scattering processes, 
since they do not change the occupation number of either 
band. Processes (a) and (b) cause the number of particles in 
the condensate band to increase, and so will be referred to as 
growth processes. The distinction between the two types of 
processes is illustrated in Fig^J Because the number of atoms 
in the non-condensate band is much larger than that in the con- 
densate band, the scattering will be dominated by processes of 
type (c), and those of type (d) will be neglected. 

These processes are described by the full quantum kinetic 
master equation obtained in QKIII, which can be used to de- 
termine rate equations for the evolution of the system. 



E. Growth Processes 



The formalism of [ |15| , |1 1||13| ] gives rise to rate equations for 
N, the number of particles in the condensate band as a whole, 
and n m , which represent the number of quasiparticles in the 
mth quasiparticle level. The derivation gives equations in the 
limit that N is sufficiently large for us to write no « N. The 
rate equations take the form: 



dn Tl 



dt 



I- ™ I growth 



'mlscatt 



(18) 



The form of h r , 



vth 



was given in QKIII, Sec.IV.E.3, and 



can be written in terms of the transition rates as follows 
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N = 2W+ I (1 
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)N + 1 



(20) 



where 



2W++(N) ((1 - e (nc(N)-n+e m )/k B T^ nm + ^ 



2\v: 



l + (N) ((1 



(21) 



,(-Hc(N)+^+e m )/k B T^ 



(22) 



the relationships between forward and backward rates can be 
shown to be 

W + {N) = e [ ^» c{N))/kT W-{N) (23) 

W 7 t + (N) = e \M-^c(N)-e m ]/k B T w -- (jy-) (24) 

W+-{N) = elc-^W+U/tB^-t^, (25) 

In these equations, the energies of the quasiparticle excita- 
tions e m are measured from the ground state energy, which is 
Mc(iV). 



1. Transitions 

There are six processes which are described by these equa- 
tions: 

i) N — ► N + 1, with no change in n. The transition probabil- 
ity for this process is W+(N) = R+(£ N ,p c (N)/h). 

ii) N — > N — 1, with no change in n. The transition proba- 
bility for this process is W~(N) = R~(£n-i, /ijv-i/fr)- 

iii) N — > N + l, withn m — * n m + l. The transition probabil- 
ity for this process is W++(N) = R+(f m , (<%+iic(N))/h). 

iv) N — > N — 1, with n m — > n m — 1. The transition 
probability for this process is W~ - (iV) = R~(f m , ( £ at-i + 
(J,N-i)/K)- 

v) N — > AT+1, withn m — ► n m — 1. The transition probability 
for this process is W+~(N) = R + (g m , (-e^ 1 + p N )/K)- 

vi) N — > N — 1, with n m — > n m + 1, The transition prob- 
ability for this process is W^ l + (N) = R~(g m , (— e r ^_ 1 + 
Pn-i)/K)- 

Here the functions f m and g m are amplitudes for the cre- 
ation and destruction of atoms in Quasipa rticle states with en- 
ergies e™, which are defined in [ p31pT| ], but whose explicit 
form will not need to be used here; £jy is the condensate wave- 
function for N atoms. The functions R^ 1 (y, to') are defined by 



R + (y,co') 



J d 3 r J d 3 K 1 d 3 K 2 d 3 K 3 d 3 k5(Aa;i 2 3(r) - u/)S(Ki + K 2 - K 3 - k)F 1 F 2 {l + F 3 )W y (r, k) (26) 
R-(y,ij')= U [ d 3 r [ d 3 K 1 d 3 K 2 d 3 K 3 d 3 k8(Au; 123 (r)-u J ')5(K 1 +K 2 -K 3 -k)(l+F 1 )(l + F 2 )F 3 W y (r,k) 



(27r) 5 n 2 



(27) 
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In these equations the following notation is used 



huj Ki (r) = 


h2K l T/ , , 

2m +VT ^ 


(28) 


AtJi23 


WKi + <^K 2 — W K 3 , 


(29) 


AK 


Ki + K 2 K 3 k, 


(30) 


u = 


47r?i 2 a 


(31) 



The function Fi = F(Ki, r) is the distribution function for 
the non-condensate particles, and the Wigner function is given 
by 



W„(r,k) 



1 



(2tt) s 



d vy* 



yir 



(32) 



The R + function is related to collisions between two non- 
condensate band particles with momenta Ki and K2, from 
which the particles leave with momenta K3 and k. The par- 
ticle with momentum k is now in the condensate band (ie: k 



is small) so that the particle with momentum K3 remains in 
the non-condensate band. The functions R ± {y, u>') therefore 
represent the rates for collisions which result in a particle en- 
tering (+) or leaving (— ) the condensate band with an energy 



2. Approximate evaluation of transition rates 

The rate factors W ± {N) contain integrals over all space 
of terms containing the product of F(K, r) terms times the 
Wigner function corresponding to the ground state wavefunc- 
tion W^ N (r, K). In practice the ground state wavefunction is 
very sharply peaked in comparison to the spatial dependence 
of the -F(K, r) functions which describe the remainder of the 
cloud of atoms. This means that in the spatial integral for the 
M /± (r) terms (from equation (p6|)) the F(K, r) terms can be 
approximated by their values at r = 0. This gives 



W + (N) 



W~(N + 1) = 



J d3K i J d " K 2 J d 3 K 3 J d 3 kS(AK)8(Au; 123 (0)- fi c (N)/h)F(K 1 ,0)F(K 2 ,0) 
x(l + F(K 3 ,0))|| Ar (k)| 2 
J d 3 K ± / d 3 K 2 J d 3 K 3 J d 3 kS{AK)S(Aco 123 (0) - Mc (iV)/fi)(l + F(K lf 0)) 
;i + F(K 2 ,0))F(K 3 ,0)|6v(k)| 2 



(33) 



(34) 



in which ^Ar(k) is the momentum-space ground-state wave- 
function, obtained from the spatial form by 



6v(k) 



1 



(2tt) 3 / 2 



d d re lkr £ N (r) 



(35) 



In [gp, progress was made by assuming that the non- 
condensate band distribution F(K, r) was given by the clas- 
sical Maxwell-Boltzmann distribution 



F(K, r) w exp 



h H /2m + Vr(r) — (j, 
kf 



(36) 



with values of T and /j, which ensure the formation of a con- 
densate once the system reaches equilibrium. Furthermore, 
in calculating the integrals in equations ( f33| ) and (Q), it was 
assumed that the range of condensate band energies was neg- 
ligible compared to that of the non-condensate band. Thus 
the range of k was negligible compared to the range of Ki 2,3 
and the integrals in Ki.2,3 were calculated over all energies 
rather than just over the non-condensate band. The function 
F(K, r) was also assumed to be negligible compared to unity. 
The result obtained was [|J 



W + {N) 



AmjakT) 2 ^ /kT 



ttTT 



(37) 



Here K\ (x) is a modified Bessel function. In almost all practi- 
cal situations the term in square brackets in the above equation 
is approximately equal to unity, and so W + is essentially in- 
dependent of N. The value of W~ (N), the rate of transitions 
out of the condensate band, can be obtained in a similar fash- 
ion to that for W + . The ratio of forward to backward rates is 
found to be given by 

W + (N) = e { ^^ c{N))/kT W- 



(38) 



which stems from the assumption of the thermal undepleted 
bath, and the definitions of W + and W~ . From this equa- 
tion and the master equation it can be seen that equilibrium is 
achieved when nc{N) = fi to order 1/N. 



3. Simple growth equation 

A rate equation for the mean number of atoms in the con- 
densate (N) (written as N for convenience for the rest of this 
section) was obtained in []5|,[lj|] 



7 



N = 2W + (N) { (l - e (Mc W-^/fc^ N+ A 



(39) 



This equation is the simple growth equation used for the sim- 
ulations of condensate growth in [|5| JTl|] . If the system starts 
with N = 0, the growth begins slowly (but at a finite rate). 
Once a significant condensate occupation is attained the term 
proportional to N becomes dominant, causing a much faster 
growth rate. The growth eventually slows as nc(N) ap- 
proaches /i and the system settles into equilibrium. This gives 
the curve an S-shape as will be shown in the following sec- 
tion. 

By using the Thomas-Fermi chemical potential ( |i"2[ ) for 
fic(N), and the Maxwell-Boltzmann form for W ± (N), the 
first simulations for the growth of a realistic condensate were 
presented in [^]. The growth equation is simple to solve nu- 
merically, for whatever number of particles is necessary (for 
example, the growth of a condensate containing 5 million 
atoms was simulated in [ph). 



4. Beyond the simple growth equation 

The derivation of the simple growth equation contains a 
number of approximations and simplifications. The major 
behavior once N becomes large should be described quite 
well by the simple growth equation, but terms which were 
neglected may have significant effects during the initial stages 
of growth. Possibly significant factors which should first be 
considered are 

i) The effect of considering all quasiparticle levels (the ex- 
cited levels in the condensate band). 

ii) The effect of scattering processes (as defined in sec- 
tion |ilD| ). 

iii) Corrections to the W /± (A r ) terms to consider the more 
realistic Bose-Einstein distribution function. 

iv) The fluctuations around the mean number. 

This paper aims to consider the effects of incorporating the 
first three of these factors into the growth equation. During the 
process of BEC formation, the spectrum of eigenvalues makes 
a transition from the unperturbed spectrum of trap levels to the 
case where the spectrum is strongly affected by the condensate 
in the ground state. The Bogoliubov spectrum of a condensed 
gas is valid in the case where the number of particles in the 
condensate, no, is so large that it is valid to write no ?« N. 
Thus, during the initial stages of condensate formation, where 
this is not true, one must use another formalism. In this pa- 
per we will consider the situation in which the interaction be- 



tween the particles is very weak, as is in practice the case. 
This means that we will be able to use the unperturbed spec- 
trum for the initial stages of condensation, and only use the 
Bogoliubov description once enough condensate has formed 
to make the effective interaction rather stronger. 

The basic formalism of [ 1 1 1 can still be carried out in this 



case, and the modification that is found is rather minor — 
essentially, we make the substitution N — > no in the chem- 
ical potential and the W + (N) , W ++ (N) functions, and set 
— > 0, since this term comes from the mixing of creation 
and annihilation operators which arises from the Bogoliubov 
method. 

By making these adjustments, and now grouping the levels 
into sub-bands of energy e m (measured now from zero, rather 
than from nc(N) as was the case for e m ), with each sub-band 
containing g m levels, the equations of motion for the growth 
processes are now 



"mlgrowth = 2W+ + (tl ) I 

"olgrowth = 2W + (n Q ) I 



£m_ti 

1 — e kT 



1 



e kT 



n + 1 



i}, (40) 
} ■ (41) 



We will make the further — possibly rather drastic — 
simplification, and neglect entirely the effect of phonon-like 
quasiparticles, which are known to comprise only a very small 
fraction of the levels normally occupied at the temperatures 
considered. Thus the excited states are now taken to be of a 
purely single-particle nature, and the condensate band is now 
described by the occupation number of the condensate level 
(the lowest energy level) no, and by the occupation numbers 
of each of the excited states n rn . In this case, the equations 
( |4C|j67| ) become the same as the equations (|i"8]-|25|), and may 
therefore be used to represent the full condensate growth pro- 



5. Evaluation of Transition Probabilities 

The value for the transition probability W + (N) found in 
the simple growth equation ( f39| ) was derived by ma king some 
rather sweeping assumptions, and as such equation ( p7| ) is re- 
ally just an order of magnitude estimate. To obtain a more 
accurate value the full Bose-Einstein distribution must be 
used for F(K, r) and the ranges of integration of the non- 
condensate functions must exclude the condensate band in 
which F(K, r) would become very large. We then have from 
d3l 



W+i - 1M)) = jirftf I Sv J d3Kl ft 3 * 2 / rf3R3 Jd 3 kS(AK)S(Auj 123 (v) - ficM/^FiK^r) 

E~>Er E~>Er E>Er E<CEr 

xF(K 2 ,r)(l + F(K 3 ,r))W 4 „ (r,k) (42) 
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with 



F(K,r) 



exp 



fr 2 K 2 /2m + V T (r)- n 



- 1 



OO 

^exp 



kT 

h 2 K 2 /2m + V T (r) - fi 



kT 



(43) 



(44) 



Again we have to make the approximations: 

i) That the spatial dependence can be neglected, so 
F(K,r) -> F(K,0) and W ?JV (r,k) -> |£v(k)| 2 . 

ii) That we can neglect the k dependence, except in 
£n, so that the only k dependence left is removed by 
/ d 3 k\£, N (k)\ 2 = 1. The integrals over Ki, K 2 and K 3 can 
then be performed, to give a final form for W + (no), found by 
Davis n of 



[log(l - z)f 



-z 2 J2lz4no)TMz,l,r + l)} 2 
r=l 



(45) 



where 



'-(no) 



(46) 



The function $ is the Lerch transcendent [h_8h, defined by 



&(x, s, a) 



OO 

E 



(47) 



This form of W + (no) gives values of about a factor of three 
greater than the previous form in equation ([37|), depending on 
the exact parameters of the system, and this gives a corre- 
spondingly faster growth than that in [||] . 

The values for W^ + (no) are more difficult to obtain. The 
W^ + (no) terms are the average of the W ++ terms for all the 
individual levels in the sub-band. The W ++ terms are given 
by similar overlap integrals as used for the W + terms, and 
for the lower energy levels in the condensate band, the over- 
lap of the wavefunction with the spatial distribution of a non- 
condensate band particle should be similar to that for the con- 
densate level. Thus it is expected that the W ++ terms should 
be of the same order of magnitude as W + (no). Progress can 
therefore be made by approximating W+ + (n ) as W + (n ). 
The effect and validity of this approximation will be investi- 
gated in the following section. 



at any time, the number of atoms in the condensate band is 
small relative to the number in the non-condensate band. The 
dominant scattering processes, and the only ones which will 
be considered, are the scattering of atoms between levels in 
the condensate band, due to interactions with non-condensate 
band atoms (see Fig.^J). These stem from terms in the full 
master equation which involve two condensate field operators 
(j>. These terms give rise to a master equation of the form, as 
shown in QKIII equation (50d): 

Plscatt = ^2 "fkmNkrn | 2 ^fcmP^| m - [X| m Xfc m ,p] + j 

mk 
'k<>m 

+ ]T j km (N km + 1) [2X km pXl n ~ [xl m X km ,p] + } 

mk 

+ E lk m M krn ^2X km pX\ m - [Xl m X km , p] + j 

km 

(48) 

which is equivalent to the master equation governing the scat- 
tering of particles by a heat bath jl9|]. 
Here the operators are defined by 



Xhm = Cl^ak 



(49) 



where a k is the destruction operator for an atom in state k with 
energy e&. As in the previous section, we treat all excitations 
as being particle-like. The rates of the processes are deter- 
mined by the factors -fkm, and the factors N km are defined 
by 



1 



exp I 



kT 



- 1 



(50) 



The last line in the master equation (|48|) represents scatter- 
ing between degenerate energy levels, which will not have any 
contribution to the time dependence of n m once the levels are 
grouped into sub-bands, and so can be ignored. 

The corresponding rate equation for n k — (a\a k ), the 
mean occupation of the fcth level, can easily be found from 
the master equation. When levels are grouped into sub-bands 
with mean energy e k , occupation n k , and with g k levels con- 
tained in the sub-bands, it becomes 

^fc|scatt 

7fem {N km n m (n k + g k ) - (N krn + l)(n m + g m )n k } 



F. Scattering Processes 



} J lik {(Ni k + l)ni(n k + g k ) - Ni k (ni + gi)n k ) 



Scattering processes in this paper also need to be included 
in the evolution of n m . Scattering between two atoms in 
the non-condensate band does not have to be explicitly con- 
sidered, since it has been dealt with in making the assump- 
tion that the non-condensate band is an equilibrated time- 
independent thermal bath. Furthermore, the scattering be- 
tween two condensate band atoms will be neglected since, 



(51) 

The transition rates j km now represent averages over the all 
the individual level transition rates which transfer an atom 
from the k sub-band to the m sub-band. The transition proba- 
bilities 7fe m can be found in a similar manner to the W + (no) 
terms, from terms of the form n\J 
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Rkm(N) = Jj" a y"rf 3 r ^^Ki ^ d 3 K 2 /d 3 k yW<f(Ki-K a -k + k^(Ki,r)(l + i?(K2,r)) 



xm(JV, r, k)W m (.N, r, k')<J(Awi2(r) - H m + Sl k ), 



(52) 



where 



n m (jv) 



Hc(N) + e m (N) 

h 

e m {N) 



(53) 
(54) 



and the rest of the notation is as was used in the previous sec- 
tion. 



1. Estimates for ■Jkm 

Explicit computation involved in calculating these factors is 
impractical, and, it will turn out, unnecessary when the scat- 
tering is sufficiently strong. We shall instead estimate these 
jkm rates by using the quantum Boltzmann approach of Hol- 
land etal. |po|l. 

By treating the excitation spectrum as given by the eigen- 
states of the trapping potential, without modification by the 
presence of the condensate, and by using the ergodic assump- 
tion Holland et al. obtained the kinetic equation (equation 
(12) in [poll) 



9fn _ \ -> c i \ 

9n — / < °e„+e m ,e p +eg9\ e mi e rai e pi Sq) 



(55) 



/,/p(l + fm)(l + fn) ~ (1 + /,)(! + f P )fmfn 



in which r = (8ma 2 Lu 2 /wfi)t. The population of a level with 
energy e n is f n , and the degeneracy of levels at that energy 
is g n . The collision kernel g(e m , e n , e p , e q ) is given from the 
overlap integrals for the states m, n, q and p. However, when 
the energies of the levels considered are spread quite far apart, 
Holland et al. found from numerical calculations that the col- 
lision kernel g(e m i a , e n ,e p ,e q ) is well approximated by the 
degeneracy <? emin . Here e m ; n is the smallest of the energies in 
the collision. In our model the energies will always be quite 
well spread, as: a) they must be in different sub-bands, each 
sub-band being quite well separated from the next in terms of 
mean energy; and b) the scattering processes we are attempt- 
ing to describe must have e m and e p in the non-condensate 
band and e q and e p in the condensate band. Thus in our model 



we may safely use g(e min , e„, e p , e q ) ps g £min . 

By summing over m and p terms (which are levels higher 
than En), me effect of all the non-condensate ('bath') levels 
on the condensate band atoms may be calculated. The kinetic 
equation for the scattering now becomes 



^2 ( l + fn)fq I J! 6 ( n > m > P > 9).9e min /p(l + fn 
q L \ "^P 

-(1 + fq)fn I ^2S{n,m;p,q)g en _ iiii f m {l + f p ) J 

\ mp / 

where the following notation has been used 

fraQyn 

6(n,m;p, q) = S en+em . ep+eq . 



(56) 



(57) 
(58) 



The terms of this equation can be simplified for the different 
possible cases. 



First Line, case e n > e q : In this case, e q 
conservation is satisfied when 



fko„ 



and energy 



(59) 



where Tiu> nq = e n — e q . The summation term in the first line 
of equation (|56[) then becomes 

g q E (i + fm)f m+ ^*g q e^- h ^ )/kT r(T) (60) 

e m >E R 

where T(T) = ^ e" 6 ™^ (61) 

e m >E R 

The approximation which has been made is that (l+/ m ) ~ 1, 
which should be acceptable since state m is of high energy (ie: 
in the non-condensate band). The calculation of T(T) requires 
knowledge of the spectrum of energies in the non-condensate 
band, which is complicated for an anisotropic trap. The form 
for an isotropic harmonic potential T(T) is easily calculated 
though, and gives 



f(T) = 



e -E R /kT 
I _ e -huj/kT 



(62) 



where u is the frequency of the potential. We will therefore 
make the approximation that T(T) r* f (T) for the determi- 
nation of 7fe m , using the geometrical mean frequency of the 
real trap as the frequency of the isotropic potential. Thus 
to = (ojxUjyUJz) 1 / 3 in equation (|62|). 
Remaining terms: Similar reasoning leads to the results 

9n ]T (1 + fp+^JU-gn J2 ^- ep)/kT (63) 

e p >E R e p >E R 

= g n e^ kT T(T) (64) 
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9q 



f m (l + f m +Hu> nq )~g q e» /kT r(T,R) (65) 



e m >E E 



9n ]T / p+ ^„(l + / p )«.g„ e ^-^")/ fcT r(T) (66) 

e p >E R 



2. Total Scattering Equation: 

The total kinetic equation governing the scattering pro- 
cesses is now given by 



2, ,2 



"mlscatt 



8ma u; 



■e^ kBT T(T) 



Y — \M9m + n m )e- hw ™«l kBT - n m (g k + n k ) 



k. < 7ii 



V — \n k (g m + n m ) - n m {g k + n k )e- huJk ^/ kBT 



(67) 



The notation k > m is now being used to mean e k > e m . 
This is the rate equation governing scattering processes, it is 
equivalent to equation (|T]) if the following transformations 
are made 



1 + N{uj nq ) -> 1 



Inq 
Inq 



8ma 2 m 2 e^ kT T{T) 

irh g n 
8ma 2 u 2 e^ kT T{T) 



nTl 



9q 



(68) 
(69) 

when n > q (70) 
when n < q (71) 



G. Rate equations including scattering and growth 

The total rate equation governing the evolution of this sys- 
tem is then given by adding equation (Eoh to equation (p7[) 



| growth ~t~ Tl m | scatt 



(72) 



and for the condensate level evolution ( |4l| ) is used in place of 
equation (f4(i|). 

It is useful here to review the major approximations that 
have been made in order to derive these equations. It was as- 
sumed that 

i) The non-condensate band is very large, so it is essentially 
undepleted, and it is in equilibrium. 

ii) The influence of collective excitations is negligibly small, 
so that the states in the condensate band are all of single- 
particle nature. 

iii) The density of states in the system is as described in sec- 
tion fic 



iv) The ergodic approximation is valid, and that states in the 
condensate band which have similar energies may be 'binned' 
together for the purpose of describing their evolution. 

v) The fluctuations of the occupation numbers around the 
mean numbers may be ignored. 

vi) The rate constants W ++ for the growth processes which 
change the occupations of the excited states in the condensate 
band are equal to the rate constant for the growth of the con- 
densate level. 

vii) The rate constant for the scattering processes in an 
anisotropic well of geometrical mean frequency lo is equal to 
that in an isotropic well with frequency lo. 

viii) Three body collisional processes may be ignored. 

III. NUMERICAL SOLUTIONS OF THE GROWTH 

The rate equations derived to describe the growth of a con- 
densate in the previous section are quite straightforward to 
solve numerically, and the solutions can be obtained in a mat- 
ter of a few seconds as opposed to other numerical solutions 
which have been very time consuming. The nature of these so- 
lutions will be discussed and comparisons will be made with 
experimental data published in [ffl]. The parameters of the 
system modelled were chosen to be the same as in the MIT 
growth experiments, that is: 

i) Using a dilute gas of 23 Na atoms, characterised by an s- 
wave scattering length of a = 2.75nm [21 3]. 

ii) With an axially symmetric ('cigar shaped') harmonic trap- 
ping potential described by the frequencies lo x = uj y = 
2ir x 82.3Hz and lo z = 2tt x 18Hz, and giving a geometrical 
mean frequency of ui = 2n x 50Hz [p]p2|]. 

iii) With a total number of atoms of the order of 10 7 . The 
numbers found in the condensate level once equilibrium is 
reached are in the range 5 x 10 5 to 1 x 10 7 , giving a con- 
densate occupation of between 5 and 30% of the total number 
of atoms [Q,^2|]. In the majority of cases the thermal bath was 
therefore depleted by only a small amount and the undepleted 
model which is used here should be a good approximation. 

iv) With temperatures in the range 0.5/^K to 



A. Results 

1. Simple Growth Equation Results 
Numerical solutions of the simple growth equation 

n = 2W + {n ) { (l - e^ c ("°)-^/ fcBT ) n + l} , (73) 

are easily obtained [||], and an example of the resulting growth 
curve is shown in Fig.^. This curve shows a characteristic 
S-shape, the slow initial growth occurs as a result of sponta- 
neous ( + 1) terms in equation ( f73| ) and then, once the oc- 
cupation becomes large enough, the stimulated growth terms 
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(those proportional to no) dominate and the growth acceler- 
ates. The condensate grows quickly, until it approaches equi- 
librium where it slows again as nc{ n o) approaches [i, giving 
the final part of the S-shape nature. 

6x 10 6 
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FIG. 6. Typical results of the simple growth equation for the 
growth of a Bose-Einstein Condensate in the MIT apparatus at a 
temperature of 900nK with \x = 43.3?lw. The initial condensate 
occupation was taken to be 100 atoms. 

The derivation of the simple growth equation contained a 
number of significant approximations and assumptions. In ef- 
fect it describes the situation in which the occupations of all 
levels higher in energy than the ground state are treated as 
time independent, giving an undepleted thermal bath in con- 
tact with only the condensate level. Furthermore, the popu- 
lations of the levels in this bath are given by the Maxwell- 
Boltzmann distribution, rather than the correct Bose-Einstein 
distribution, and no mean-field effects are introduced. The 
simple growth equation is therefore just a first approxima- 
tion for the growth, and was merely intended to give an order 
of magnitude description of the growth process. Indeed, as 
shall be shown a subsequent section, experimental measure- 
ments [Q] have shown that, although it gives the correct order 
of magnitude for timescale of the growth, the simple growth 
equation does not describe the rate of growth to any closer 
than a factor of about 3. 



2. New Model of Growth Processes 

In order to improve the description of the growth over that 
of the simple growth equation, the first improvement which 
will be made is the more accurate calculation of W + , us- 
ing the full Bose-Einstein distribution. However, for lower 
energy levels the equilibrium populations determined by the 
Bose-Einstein distribution are very large. Having such large 
populations in these levels will obviously not be a good model 
of a system rapidly cooled from a point where fj, was negative 
to a region where /i has become positive, since the changes 
in populations of these levels required during this process are 
so substantial that a fairly long time will be required for the 



levels to come to equilibrium. It is therefore unphysical to 
consider a situation where the low energy levels have reached 
their equilibrium populations before the condensate level has 
even started to evolve. For these reasons, it is not consistent 
to simply use the Bose-Einstein distribution to find W + in the 
model considered by the simple growth equation. 

In order to develop a consistent description there need to be 
a number of lower energy levels with time-dependent popula- 
tions considered, forming the condensate band. Furthermore, 
the energies of these levels will be increased by the growth 
of the condensate, due to increased mean field interactions. 
To describe this situation the model described in the previous 
section must be implemented. Considering at first on ly the 
'growth' processes, the evolution was found in section HE 4 
to be given by 



| growth 



n \ 



vth 



2W+(n ) { " 



1 - e i*t 



1 — e kT 



n m + flm } , (74) 
no + l}. (75) 



where, for reasons given in the previous section, we approx- 
imate W^ + (no) by W + (n^). This form of evolution is es- 
sentially the simple growth equation applied to several energy 
levels. 

It is now possible to perform a more accurate calculation of 
W + , using the Bose-Einstein distribution to describe the pop- 
ulation of levels above Er (ie: in the non-condensate band), 
and summing only over the non-condensate band levels. The 
new form of W + is given by equation (pBj). 

Sample solutions to the coupled differential equations ([74j) 
and ( [75| ) are shown in Fig .[7]. The rate of growth of the conden- 
sate has increased substantially (generally by at least a factor 
of 3) over that predicted by the simple growth equation, this 
is due to the more accurate calculation of W + (no). However, 
the shape of the growth is still essentially the same as that 
given by the simple growth equation. 

It can be seen from Fig.^ that the lower energy levels also 
experience very substantial growth in this model. Indeed the 
occupations of some of these levels can exceed the condensate 
occupation substantially before relaxing back to their equilib- 
rium values. This is of course not a realistic scenario, and is 
certainly not one that has been observed in any experiments. 
Note that the number of sub-bands used in this figure (and 
most of the other figures presented) is substantially fewer than 
would normally be used. Generally the number of sub-bands 
required is about 20-50 depending on the exact parameters, 
however these cannot be well distinguished from each other 
in a graph. Therefore, throughout this paper, most of the de- 
pictions of the growth of sub-bands will show only a few of 
them. 
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FIG. 7. Typical results of the new growth equations for the growth 
of a Bose-Einstein Condensate using the same parameters as in Figjfj. 
The condensate level (bold) reaches the equilibrium population of 
5 x 10 6 atoms, the other lines represent the evolution of the popula- 
tions of the other sub-bands in the condensate band. Note the much 
faster growth than in Fig.^j. 

Once scattering processes are included (see below), reduc- 
ing the number of sub-bands used causes the growth to be- 
come slower. It also causes the model to become less realistic, 
since an individual level may then be described by an average 
energy quite different from its actual energy. If the number of 
sub-bands is increased, an asymptotic limit to the speed of the 
growth is reached, however this is also unrealistic since now 
some sub-bands contain only fractions of individual energy 
levels. The choice used in practice is such that in their fi- 
nal equilibrium states the sub-bands have widths of Tito. This 
choice is close to the asymptotic limit, and ensures at least 
three individual levels are contained in the first sub-band. 



3. Inclusion of Scattering Processes 



If we include scattering processes, as given by (^), the pic- 
ture is dramatically changed. Solutions for the resulting evo- 
lution equations (|72|) are shown in Fig.|| This figure shows 
that the scattering has two main effects. Firstly, the initiation 
of the condensate level growth occurs much more sharply, this 
gives a substantial change to the shape of the growth, which 
has now lost much of the S-shape nature that previous solu- 
tions had shown. The speed of the growth after the initiation 
is changed little by the inclusion of the scattering processes, 
since in this region growth is completely dominated by the 
growth processes. The second effect is that the populations of 
the excited states no longer exceed plausible levels. 



6x10 



4x10 



2x10 



2x 10 5 



S 

Z 



% lx 10 5 




120 160 200 
Time (ms) 

FIG. 8. Typical results of the total growth equation (including 
scattering processes) for the growth of a Bose-Einstein condensate 
using the same parameters as in Fig^. a) The growth of the con- 
densate (ground state) occupation. The inset shows the same growth 
predicted by the total new model (solid) compared to the growth pre- 
dicted without considering scattering processes (dashed), and that of 
the simple growth equation with an initial population of 100 atoms 
(dotted), b) The growth of the explicitly considered sub-bands. The 
condensate level itself is the bold curve. The top five sub-bands are 
those in the non-condensate band which were explicitly considered. 
Note the different timescales on the graphs. 



The reasons behind these changes are interlinked. With- 
out scattering, all of the levels in the condensate band start to 
grow, and at quite similar rates. The difference in the growth 
rate between the very lower energy states and that of the con- 
densate level is particularly small. Because of the degeneracy 
of states in the lower sub-bands a large population can form 
in them, which can become very much larger than the con- 
densate level population. Once any one sub-band acquires a 
sufficient population the stimulated term in the growth process 
begins to dominate and the population increases even further. 
In the absence of scattering processes, the only way in which 
the excess population in these states can be transferred to the 
condensate level, where it will be found in the equilibrium sit- 
uation, is by a transfer back to the non-condensate band fol- 
lowed by another collision which transfers it directly to the 
condensate level. 

If scattering processes are considered, atoms may now be 
transferred directly between different levels in the condensate 
band in a collision. Any excess population in the excited states 
can then be quickly transferred out of the state before the stim- 
ulated growth process becomes too dominant. 
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With the inclusion of the scattering processes, the effects 
of two important approximations must be considered. In the 
derivation of the equations governing the scattering processes 
it was assumed that the value of T(T) = ^ e >Er eT e ™l kT 
was equivalent to that for an isotropic harmonic oscillator 
f (T), which affects the scattering rates. 

The second important approximation was made in the solu- 
tion of the equations governing the growth processes, where 
the W+ + (no) terms were assumed to be equal to W + (n ) 
which has a value given by equation (45). This approxima- 
tion has no effect on the growth rate of the condensate level if 
the scattering terms are not considered, since then there are no 
interactions between different levels in the condensate band. 
The effects of these two approximations are discussed below. 
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FIG. 9. The dependence on the rate for scattering processes of 
the condensate level growth. Growth curves from left to right had 
factors T(T) of lOOf (T), f (T), f (T)/2, f (T)/10, f (T)/100 and 
zero respectively, where T(T) is the factor for the isotropic harmonic 
potential given by equation (M). These results were obtained using 
the same parameters as used in Figj8| 



4. Effect of the Scattering Rate Approximation 

The magnitude of the scattering rate was assumed to be 
equal to that for an isotropic harmonic oscillator in the deriva- 
tion of equation (^). In Fig|| the growth of the conden- 
sate level is shown for several different scattering rates. The 
growth slows slightly if the rate is decreased. If the rate is in- 
creased, the growth becomes faster until it reaches an asymp- 
totic limit at which point it is the rate of the growth processes 
which determines the speed of growth. The results show that 
the overall growth changes by only a relatively small amount 
(and certainly smaller than present experimental uncertain- 
ties in growth experiments) provided that the rate is within 
two orders of magnitude of that for the isotropic trap. Since 
it seems unlikely that the corrections due to the anisotropy 
would change the rate by much more than one order of mag- 
nitude, the approximation of using the isotropic rate factor 
seems to be valid. It is interesting to note that a scattering rate 
of only about 1% of the isotropic case is usually sufficient to 
prevent the occupations of the sub-bands becoming very large 



as they do in the absence of any scattering 
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FIG. 10. The effect on the rate of growth of the condensate when 
the values of W^ + (no) were changed. From left to right the curves 
correspond to values for W+ + (no) of 10W + (no), 2W + (no), 
W + {n ), W + {n )/2 and W + (n )/10. Conditions correspond to 
those in Fig.0. 



5. Effect of factors 

The effect of the approximation W+ + (no) ss W + (no) is 
shown in Figjl(] which shows the behaviour of the growth for 
different values of the W+ + (no) terms. Changing W+ + (no) 
changes the rate of growth of the excited states, and thus 
changes the probability of atoms being scattered from the ex- 
cited states into the condensate level, giving a corresponding 
change in the overall growth of the condensate. The results 
show that provided W+ + (n ) lies in the range 2W + (n ) > 
W+ + (n ) > W + (n a )/2 then the growth rate does not 
change significantly (compared to experimental uncertainties 
and the change caused by using an accurate scattering rate). 
However, outside of this range the growth is altered consid- 
erably. The expectation is that W^ + (no) will lie in the de- 
sired range, since it is an average over quantities similar to 
W + (no), and as such should be of the same magnitude. 

The 'standard' approximations for the rate constants that 
will be used in the following results will be W^ + (no) — 
W + (n ) andr(T) = f(T). 



6. Initial Conditions 



A problem which has to be considered in solving the total 
rate equations ( |72| ) is the determination of the correct initial 
conditions. Because the non-condensate band is assumed to 
be a thermal bath of atoms at equilibrium, the initial popula- 
tions for the explicitly considered levels in this band are found 
from [cxp((i? — n)/kT) — Obviously this cannot be used 
to give the initial populations in the condensate band when the 
model is attempting to describe growth of the condensate. It 
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is not immediately obvious what the appropriate initial condi- 
tions should be for the condensate band. In Fig. H] four differ- 
ent initial conditions are shown 

a) no initial population in the condensate band. This is the 
most artificial of the four possibilities presented. 

b) initial population of zero in the condensate level and ex- 
cited state occupations given by a linear dependence on energy 
rising to match that of the non-condensate band at Ef>. 

c) initial populations of condensate band states all equal to 
the value of [exp((E R - n)/kT) - l]" 1 . 

d) initial populations given by a linear extrapolation of the 
non-condensate populations, meeting [cxp((E — /j,)/kT) — 
l]" 1 at Er tangentially. 

In Fig.[j~2|the growth curves corresponding to each of the dif- 
fering initial conditions are shown. The different initial con- 
ditions can be seen to have little effect on the shape of the 
growth, the main effect is really just a small shift in the ini- 
tiation time. This effect is generally quite small compared to 
the effects of changing T\T) and W+ + (no) as seen in the 
preceding sections. The fact that the initial conditions can be 
changed by so much and yet have little effect on the growth 
curve is due mainly to the inclusion of the scattering terms. 
The scattering terms very quickly cause the levels in the con- 
densate band to come to a kind of quasi-equilibrium from 
whatever initial state they are put in. Thus we conclude that 
an exact knowledge of the intial conditions is not important. 

Since the exact initial conditions do not seem to be im- 
portant, in all further calculations initial condition (c) will be 
used, as it is about mid-way between the extremes of cases (a) 
and (d). 
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FIG. 12. Growth curves determined using the initial conditions 
described in the text. The curves, from slowest to fastest growth, 
were given by initial conditions a-d respectively. Conditions are the 
same as for Figj8| 



7. Modifications of the Energies of the Sub-bands 

The previous results were all obtained by taking the mean 
field effects due to the condensate into account in the manner 
described in the previous section. That is to evenly distribute 
the energies of the sub-bands between the fixed upper limit of 
Ea and the no-dependent lower limit /J,c{ n o)- Thus all levels 
below Ea (fixed at 2nc{ n o,f) where uqj is the final occu- 
pation of the condensate level) are modified at all times. This 
artificial model can be improved on, since the number of lev- 
els affected by the condensate depends upon its occupation. In 
other words, when there is only a small condensate present it 
has a significant effect on only the lower sub-bands, while the 
upper sub-bands are essentially unmodified. An improved (al- 
though still a little artificial) method is the following: all lev- 
els below 7"c( ri o) ( wnere 7 is an arbitrary parameter) will be 
compressed to fit between this upper limit and the lower limit 
of nc(n>o)', the energies of all levels higher than 7/xc(no) will 
not be altered. This is similar to the previous model, with the 
alteration that the upper limit of the levels whose energies are 
changed is no longer fixed, but instead rises with increasing 
no. The energy of a level e m is then given by 



Mc(«o) + e£ 



1-i 



fore°„ > 7u c (n ) 
fore m < 7Mc(™o) 



(76) 



FIG. 11. Depiction of the four different initial populations trialed, 
as described in the text. 



where e® n are the unchanged energy levels, as given by the 
eigenstates for a non-interacting gas in a harmonic poten- 
tial. The energy spectrum of the sub-bands given by this new 
model is shown in Fig|l3|, and compared to the previous en- 
ergy spectrum. In the older model used so far in this paper 
all energy levels in the condensate band are changed by the 
growth of the condensate. However, the extent of the con- 
densate band is determined by the final occupation number 
of the condensate (i.e., by the final value of /xc(^o))- In the 
new model, the levels are modified in a more consistent fash- 
ion, the energy levels for any given condensate population are 
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determined by the condensate population at that time, which 
appears to be a more logical approach. 



In section [IC2 it was shown that a fair estimate for Er 
was a value of 2//c(^o,/)- The value for 7 will normally be 
taken to be 2, so that at equilibrium all energy levels below 
Er will be modified. 
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FIG. 13. The energy spectrum of the sub-bands as a function of 
no, modified to account for the mean field interactions with the grow- 
ing condensate. In (a) the levels are modified by changing all levels 
in the condensate band, as has been used to obtain the previous re- 
sults. In (b) these effects are accounted for by the new model, ex- 
plained in the text. The lowest energy level is the condensate level, 
whose energy is fj,c(rio). The dotted line represents 2pc{no), the 
maximum energy at which mean field corrections are assumed to 
be significant. This figure shows the evolution of the energies for a 
temperature of 900nK and with a final (equilibrium) condensate oc- 
cupation of 5 x 10 5 . The few levels whose energies never change 
are those levels in the non-condensate band (above En) which are 
explicitly considered. 

The effect of this new model, compared to the earlier ap- 
proach, on the overall growth is very small, smaller in fact 
than any of the other effects discussed in this section, and it is 
barely discernible. This would seem to indicate that it is the 
mean-field effects on the energies of the lower few sub-bands 
which are important when considering the growth of the con- 
densate, since in both the old and new models the lower levels 
experience quite similar changes. Since the mean-field effects 
of the condensate on the higher levels seems to have only a 
very small effect, the precise value of Er would appear to 
have little effect so long as it is reasonably high. The value of 
2/j,c(no.f) seems to be a good value, since it is high enough 
that the energy perturbations of higher energies do not have a 
large effect on the growth, and it is low enough that the major- 
ity of atoms are found at higher energies, giving an undepleted 
thermal bath. 



B. Comparison of Results with Experiment 

At the present point in time there has only been one pub- 
lished experiment which has investigated in any detail the 
growth of a Bose-Einstein condensate. The experiment was 
performed at MIT using a trapped gas of 23 Na atoms, and the 
results were published in 1998 in ref fl]. 



1. MIT Experimental Method 

The MIT experiments were performed in the following 
way: The gas of atoms was confined in an approximately har- 
monic magnetic trap. It was then cooled using laser cooling 
and evaporative cooling techniques to a temperature slightly 
higher than the critical temperature necessary for the forma- 
tion of BEC. At this point the system is essentially in ther- 
modynamic equilibrium, and then it is suddenly put into a 
non-equilibrium configuration of lower energy, by means of 
a rapid evaporative cooling 'cut' which removes all atoms in 
states above a certain energy in a time of about 10ms. The sys- 
tem is then left to relax to equilibrium with no further cooling. 
The cut will have brought the temperature of the gas below the 
critical temperature and so, to reach equilibrium, a conden- 
sate will form. The formation of the condensate is observed 
at several stages during the evolution by the means of phase- 
contrast microscopy. 

This method attempts to achieve, probably as closely as is 
realistically possible, an almost thermalised bath in contact 
with a condensate band, as has been assumed in our theo- 
retical treatment. The cut which removes the higher energy 
atoms causes the wings of the energy distribution to be trun- 
cated. Experience with solutions of the quantum Boltzmann 
equation and related methods [ p8|j2Q| , ^0| ] shows that the effect 
of this will firstly be felt by the higher energy atoms remain- 
ing. The higher levels will therefore be expected to thermalise 
more quickly, with thermalisation gradually moving through 
to the lower energies. Thus at some point after the cut the ma- 
jority of the atoms will be approximately in equilibrium, with 
the lower energy atoms still in very non-equilibrium states. 



2. MIT Experimental Results and Theoretical Comparisons 

Phase-contrast microscopy produces two dimensional im- 
ages of the system with an intensity proportional to the col- 
umn density of the system. From data of this type tempera- 
tures, total numbers of atoms, and condensate level occupa- 
tions were extracted by the MIT group and were presented 
in[|. 

The results found in [|4| were that condensate growth took 
on the order of 100-200ms depending on the exact conditions. 
It was found that the growth could be well fitted by a solution 
of the simple growth equation (J73[), which can be put in the 
approximate (but in practice very accurate) form 



h Q = «ino[l - (n /n ,f) 2/5 ], 



(77) 



where, again, uqj is the equilibrium condensate population. 
The solutions to this equation exhibit the S-shaped growth 
profile of Fig.||. 

The conclusion drawn in Q was that a curve of this shape 
was evidence for the importance of bosonic stimulation in the 
growth processes, since a purely relaxational process would 
be described by solutions of 



h = K 2 (noj - no). 



(78) 
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However, the rate constants Ki found by the MIT group by 
fitting to the data obtained did not agree to better than an or- 
der of magnitude with the predictions of the simple growth 
equation of 



Ki = 2W + {n a ) 



JL 
kT' 



(79) 



The simple growth equation did well to give the correct or- 
der of magnitude rates, however most of the predicted rates 
for higher temperatures seemed to be too small by about fac- 
tor of three. As the temperature decreased the discrepancy 
increased, the predicted rates became slower, whilst the ex- 
perimentally fitted rate constants became larger. 
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FIG. 14. Comparison of theoretical growth curves with data ex- 
perimentally obtained at MIT (dots). The measured temperature for 
the experimental data was 1200nK. In the main part of the figure the 
solid line shows the theoretical fits with standard rate constant ap- 
proximations at 1200nK, the dashed line shows the theoretical curve 
obtained by using f (T) = r(T)/10 and W+ + (n ) = W + (n )/2. 
The dotted line shows the growth predicted by the simple growth 
equation of [Q] (equation (|3S|)), which is shown on a larger scale in 
the inset. The initial condensate populations for each curve were 
set to 5 x 10 , the experimental value at t = 0. The origin in the 
time axis represents the time at which the quick cooling 'cut' in the 
experiment was finished. 



In thegrowth experiments statistical uncertainties are esti- 
mated [ p9| ] to be 10% for relative number measurements, 15% 
for temperatures. Systematic uncertainties are estimated as 
20% in absolute number measurements, and 8% for tempera- 
ture measurements. Condensate occupations of less than 10 5 
atoms could not be discerned against the background of the 
thermal vapour cloud [f|] . 

As Fig.[l4] shows, the growth predicted by our model is 
quite a good fit to the data, and the order of magnitude is 
certainly predicted well. This is a substantial improvement 
to the growth predicted by the simple growth equation at the 
same conditions which gives growth over about 1.5 seconds, 
as opposed to the experimental results of about 0.15s. At the 
stated parameters the theoretical fit could still be improved. 
The dashed line in Fig.[l4| indicates that by adjusting the scat- 
tering rate and values of (as was discussed in previous 
sections) a better fit may be obtained. 
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FIG. 15. Comparison of theoretical growth with MIT experi- 
mental data. The experimental data is the same as that used in 
FigJI^. The theoretical curves were determined by using tem- 
peratures of lOOOnK (solid) and 850nK (dashed). The rate con- 
stants were taken to be their standard values of T(T) = T(T) and 
WX + { n o) = W + (no). Again, the initial condensate populations 
for each curve were set to 5 x 10 4 . 




In contrast, the solutions to the growth model presented in 
this paper no longer show the S-shape, but are in fact closer 
in shape to the solutions of equation ([78]). These curves were 
found by MIT to describe the data quite poorly if the growth 
started at time t = (the time at which the cut finished). How- 
ever, if an initiation time was allowed before the growth be- 
gan, such solutions became quite close fits, although they still 
did not describe the initially slow growth giving the S-shape to 
the growth. The results obtained in this paper show that, while 
the general shape is that of solutions to equation ([78]), there is 
also an initiation time present before the growth begins. 

As a specific comparison with experiment, Fig.[l4] shows 
the comparison between the growth curve predicted here and 
the experimental data. The data is similar to that in [Q], but 
was not actually published. It was provided by MIT as being 
the growth of a condensate with an equilibrium population of 
about 1 x 10 6 atoms, at a temperature of 1200nK [E2fl. 



The temperature plays a sensitive role in these comparisons, 
since the rate of growth is quite sensitive to temperature. Fur- 
thermore, as will be seen in the following section, the fitting 
method may play a significant role. In Fig.|l5] the theoretical 
curve is plotted (using the standard values of the scattering 
rate and W^ + ) for two lower temperatures. It can be seen 
that the fit is very good for the 850nK results, and not quite so 
good, but still quite close, in the case of the lOOOnK curve. 

This highlights a difficulty in comparing theoretical predic- 
tions with experimental measurements of condensate growth. 
The spatial density distribution of the thermal cloud (from 
which the temperature is measured) changes only slightly with 
temperature, whereas the growth rate is quite strongly depen- 
dent on the temperature. There are other problems as well, the 
most prominent of these being the difficulty of experimentally 
determining the condensate occupation. The spatial distribu- 
tions of particles in the first few excited states are quite narrow 
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functions, and they can overlap the condensate level distribu- 
tion significantly. It then becomes extremely difficult to dis- 
tinguish condensate level atoms from low energy excited state 
atoms in experimental measurement. If the measured 'con- 
densate occupation' was in fact the occupation of the lowest 5 
levels (for example) this would alter the growth curve in some 
important regions. The main effect would be that the growth 
would appear more gradual during the early times, and would 
not seem to have such a sharp initiation. This is a possible 
explanation for the S-shape that was found in the MIT data. 
Clearly what is needed are theoretical predictions for the over- 
all spatial density distributions during the growth, rather than 
merely occupation numbers for the various states. This will 
be addressed in the following section. 
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FIG. 16. Comparison of theoretical growth (thin lines) with 
curves fitted by MIT to experimental data in |Q| (broad lines). The 
MIT curve fits are represented as broad lines to indicate that they 
are fitted curves with unknown (but probably substantial) uncer- 
tainties, a) Theory T = 830nK, no,/ = 7.6 x 10 6 ; Experi- 
ment T = 810 - 890nK, n 0J = 7.5 - 7.85 x 10 6 ; b) Theory 
T = 590nK, n j = 2.3 x 10 6 ; Experiment T = 580 — 610nK, 
no,/ = 2.1 — 2.5 x 10 6 . The initial populations were treated as free 
parameters to best match the initiation of the growth with the MIT 
curves. 

The majority of the data contained in was presented in 
the form of rate constants for fits to experimental data of the 
type described by equation (|77|). Fig.|l6|shows the comparison 
of our theoretical growth curves as compared to solutions of 
equation ( |77| ) using a selection of the MIT fitted constants n\. 
The figure shows good overall agreement with the experimen- 
tal data. The agreement is better at higher temperatures. At 
the lower temperatures reached in the experiment the agree- 
ment is less good. 

However, the results of this paper still show an overall de- 
crease in the the rate of growth with decreasing temperature. 
This is the opposite trend to that experimentally observed. 
This could partly be due to uncertainties in extracting num- 
bers from the experimental data. Another possible explanan- 
tion is that in order to cool to such low temperatures from just 
above the critical temperature in 10ms a large proportion of 
the atoms must be removed. This will give rise to a system 
in a highly non-equilibrium state, the relaxation from which 
may be inadequately described by our model. 



IV. SPATIAL DENSITY DISTRIBUTION OF A 
CONDENSATE SYSTEM 

In the MIT experiments into condensate growth [^|] the data 
is collected by phase-contrast microscopy [ p4p5| ], the mea- 
surement of change in the phase of light after it has passed 
through the vapour cloud. The result is a two dimensional plot 
of the column density integrated along the third dimension. As 
discussed in the previous section, the extraction of the popu- 
lation per energy level from this experimental data is compli- 
cated. Therefore, in order to more easily enable comparison 
with experiment, it is desirable to obtain from our theoretical 
results predictions of the spatial distribution of the condensate 
as it grows. When the system is in equilibrium this distribution 
is well known, but this is not the case during the condensate 
growth that we are interested in. 



A. Semiclassical Phase-space Description 

In order to convert the results of the model into spatial dis- 
tributions, the spatial probability distributions for each energy 
band need to be found. The distribution for the ground state 
Pg S (r) is well described by the Thomas-Fermi approximation 
for the wave function: 

P gs (r) = \iPtf(v)\ 2 

= —T- \Mo(no) - V T (v)] 9 (/i C (n ) - V T (r)) , (80) 
Ann a 

where Vr(r) is the trapping potential, and 6{x) is the step 
function. This is a very good approximation to the shape of 
the condensate when no is large, failing only at the very edge 
of the condensate where the numerical solution to the Gross- 
Pitaevskii equation vanishes smoothly. 

A description of the spatial distributions for each of the 
higher energy sub-bands can be found by using a semiclas- 
sical phase-space approach. The cumulative number of states 
below an energy E is given by 



(81) 



where V^cM is the effective potential experienced by an 
atom. This gives a density of states 



9(E) = 



dN{E) 



dE 



dP S E 



p2 

2m 



(82) 



Using the local density approximation, for an energy sub-band 
with energy e m and width Ae m , the average spatial distribu- 
tion atoms in the band (averaged over all the wavefunctions 
of all the states in the band) may be obtained by removing the 
d 3 r integral and integrating over <i 3 P. This gives 
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FIG. 17. Sample radial spatial density distribution for condensate 
having reached thermal equilibrium, a) Solid — total density distri- 
bution, dashed — density distribution due to excited states only, b) 
Density times radius squared. These results are for the trap parame- 
ters used in the MIT growth experiments | 



1. Resulting Spatial Density Distribution 

We will use the semi-classical distributions for the excited 
states, and the Thomas-Fermi wavefunction for the conden- 
sate. The total radial probability distribution may be calcu- 
lated by normalising each level, and then summing over all 
the levels. The distributions are normalised so as to give 
the appropriate population in each level. For the bands be- 
low En these populations are given by the numerical results 
of our model. For those above Er the population is deter- 
mined by the Bose-Einstein distribution function F(E) = 
[cxp((_E— n) /kT) — This gives the spatial density distri- 
bution for the whole condensate system in three dimensions. 

An example of the resulting spatial density distribution for 
a system at equilibrium is shown as Fig.|l7[ In the first part of 
this figure the total density is plotted as a function of radius, as 
well as the density due to the excited states only. The density 
due to the excited states can be seen to be significantly de- 
creased in the region of the condensate. In Fig.n7p the density 
multiplied by the radius squared is plotted, this is proportional 
to the total number of atoms found at any radius (due to the 
three dimensional nature of the distribution). From this it can 
be seen that, even though the density in the centre due to the 
condensate is much larger than elsewhere, the majority of the 
atoms are still found in the surrounding vapour cloud (as was 
assumed in the derivation of the model). 

When the system is in equilibrium these results can be 
checked. At equilibrium the non-condensate spatial distribu- 



tion can be obtained directly from equation (p2|). The total 
number of non-condensed atoms, M is 



M= I g{e)- 

'0 



1 



-de 



e {e-n)/kT _ I 

which corresponds to the local density form of p(r) 



P (r) = ^[ d e[ d *r 6 ( E - P2/2m - V ^ 



(85) 



d 3 P < exp 



; ( £ -/x)/feT _ i 



£ + ^S(r)-M 
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- 1 



(86) 
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Carrying out the integral over momentum space gives 



27rftV V kT 

oo 

where G a (z) = 



(89) 



A comparison of the non-condensate density obtained using 
equation ( p8[ ) and that calculated by summing equation ( p4[ ) 
over 1200 levels is shown in Fig|l8|. It can be seen that the 
agreement between the two methods, for this equilibrium situ- 
ation, is good. The agreement improves if more energy levels 
are included in the sum over equation (|84j), and using about 
2000 levels gives very good agreement. Our semiclassical 
method therefore shows good agreement with the expected 
distribution at equilibrium, and this is the only case in which 
we can be certain of the theoretically correct result. 

In obtaining the results in this section we have only con- 
sidered effects of the mean-field repulsion due to the con- 
densate atoms, both on the thermal vapour cloud and on the 
condensate itself. In order to be truly consistent the model 
should also include the effects of the mean-field repulsion 
due to the thermal cloud on the system. This is quite easily 
achieved mathematically, but it does increase the computa- 
tional time by a very substantial amount. However, recently 
Naraschewski and Stamper-Kurn [^6j] compared the density 
distributions obtained both with and without considering the 
mean-field repulsions of the thermal cloud, for an equilibrium 
condensate system. They found that the overall density dis- 
tributions for the two cases were practically indistinguishable 
except for a very small deviation at the edge of the conden- 
sate. Furthermore, Holzmann, Krauth and Naraschewski fiFfy 
found that the semiclassical density distribution (including the 
mean-field repulsions by the thermal cloud) gives excellent 
agreement with exact Quantum Monte Carlo simulations for 
dilute gas condensates in equilibrium. Thus we expect that 
the density given by our semiclassical method is a very good 
description of the realistic system. 



19 



a 

-3 



-a 
a 

o 

o 

a 
o 
Z 



a) 


i 




II 


V 


- 




V 






V 






V 






1 [_ 





50 



100 



150 200 
Radius, r (um) 



250 



300 




100 150 200 
Radius, r (urn) 



300 



FIG. 18. Comparison of the non-condensate density distribution 
at equilibrium calculated a) from the sum of 1200 individual energy 
level distributions (dashed), and using the G3/2(z) function (solid); 
b) as for (a) but using the sum of 2000 energy levels. 



2. Description of a Realistic Experimental System 

The next step that must be performed, in order to compare 
with the experimental results, is to perform a column integral 
along one dimension, yielding a function of two spatial di- 
mensions only. This is relatively straightforward to achieve 
numerically. 

Finally the asymmetry of the real traps must be taken into 
account. The previous arguments for the exact non-interacting 
wavefunctions assumed the use of a spherically symmetric 
trap, whereas the realistic traps used are strongly anisotropic 
in the z dimension. However, in all the previous semiclassical 
arguments the only effect is to change Vr from mw 2 r 2 /2 to 
mui 2 y (x 2 + y 2 )/2 + muj 2 z 2 /2 where now 
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U)/X 

UJ = (ui x L) y U) z ) 1 / 3 ' 

A = (oj x /oj z ) 2/3 . 
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original form of the potential by scaling the dimensions used. 
Defining 



z = z/X 
x = xVx 

y = y^, 



(94) 
(95) 
(96) 



the potential now returns to the form Vr(r) = mw 2 r 2 /2, 
where the scaled radius is given by f = ^/x 2 + y 2 + z 2 . 
Equation (|4|) can now still be used, and the column density 
integration performed as for an isotropic potential, and the re- 
sulting function of x and z needs only to be rescaled to recover 
the answer in terms of x and z. This scaling also affects the 
numerical column integration, with the effect that the result 
needs to be scaled overall by a factor of A~ 3 / 2 (assuming the 
integration is in one of the two shorter dimensions and that the 
result is in the asymmetric (longer) dimension, as is the case 
for the 'cigar' geometry traps of MIT). 



B. MIT Fitting Method 

In the MIT experiments, the data for the numbers of con- 
densate atoms and temperatures are obtained from fitting to 
the density profile obtained. The raw data obtained is in the 
form of two dimensional images and, although these can be 
fitted, in order to save time and computational resources the fit 
was mostly performed only to a one dimensional slice through 
the center of the condensate fizSj. This density profile is fitted 
from a function formed by the combination of a condensate 
density profile and non-condensate profile. The MIT fitting 
procedure (presented in [[0]]) is a phenomenological proce- 
dure, only loosely connected with fundamental theory. It is 
based on two observations: 

i) The behavior of the thermal cloud in the wings of 
the profile is almost independent of the chemical 
potential — its behavior in the wings can thus be used 
to determine the temperature of the vapor. 

ii) The center of the profile is dominated by the conden- 
sate, and this — after subtracting the contribution of the 
vapor — can be fitted to the parabolic Thomas-Fermi 
profile, to determine the condensate chemical potential 
and thus the number in the condensate. 

In order to execute the second procedure, some estimate of 
the vapor contribution has to be given in the region where the 
condensate is present. This is done by fitting the profile to the 
sum of a Thomas-Fermi parabola and a zero chemical poten- 
tial Bose-Einstein distribution. 

Thus, using our methodology, this would be as follows. Set- 
ting the chemical potential to zero and integrating along one 
dimension gives (for a slice where x — 0) 



Pint (2, a; = 0) = / p(z,y,x = 0)dy 
J — 00 

m(kT) 2 



Due to the harmonic nature of the trap, one can recover the 
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This function can then be used to fit the wings of the distribu- 
tion to obtain T. 

The column density function due to the condensate, pi n t,c, 
is then obtained by integrating the Thomas-Fermi wavefunc- 
tion given by equation (pfl) over one dimension 



Pint,c{z,X = 0) 
Vi( z ) 

■Vi( z ) 



Hc(n ) m 2 2 , 2 2n 

— Yu^ z +UJ y y) 



dy 



where yi(z) 
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and where u = 4:TTh 2 a/m, giving 
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The dens ity profile can therefore be fitted by the sum of equa- 
tion (101) and equation I 



In the measurement procedure used, the phase shift data 
was not calibrated independently. Instead, the raw data were 
fitted using semiphenomenological fit functions and an arbi- 
trary scale factor to related the measured phase shif to the col- 
umn density. 

However, the actual procedure used is somewhat more gen- 
eral, since the MIT procedure uses the trap constants as further 
fitting parameters, and these are allowed to be independently 
chosen for the condensate and the thermal cloud. This is done 
since finite resolution effects of the imaging system may have 
a rather different effect on the very narrow condensate pro- 
file from that on the broader vapor profile. In practice, for 
similar reasons the scale factor relating column density to the 
phase-contrast signal is allowed to be different for vapor and 
condensate. The number of condensate atoms no can be ob- 



tained from fitting directly to equation (101), or, by a simpler 



method used for most cases, determining the extent of this 
condensate distribution. In most cases, since it is computa- 
tionally simpler, number of condensate atoms was obtained 
from the spatial width of the condensate zq. These are related 
by 



Pint.cOzO) x = 0) = when z% = 
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(102) 



thus no oc Zq. 

In our fits we have assumed the scale factor to be the same 
for condensate and vapor, since our model is physically based, 
rather than being semiphenomenological. Hence, in our pro- 
cedure, there are three free parameters to be fitted: T, no and 
a scale factor. 



C. Comparison of Fitting Methods 



In Fig. [19] we show a least squares fit of a zero chemi- 
cal potential plus a Thomas-Fermi condensate profile for an 
equilibrium distribution. The data used was obtained from 
Fig.2 of [Q], for which the MIT group obtained values of 
T = 8OO11K with 9 million condensate atoms present in equi- 
librium after 160ms jj22|]. Also shown is the contribution to 
the distribution due to the non-condensate atoms only. Note 
that the z axis has been scaled by a factor of 3.2 over that 
in [[(J due to a typographical error in the paper p2|]. 

In this fit, as for all the following cases, the density scale 
factor (necessary to convert from the arbitrary unit scale re- 
sulting from the experimental measurements) was determined 
by the least squares fitting procedure. Since we use the T and 
no determined by the MIT group using their procedure, this 
scale factor is the only free parameter in this fit. 
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FIG. 19. Equilibrium data obtained from MIT [Q] fitted using the 
MIT procedure as described in the text (solid smooth curve) and 
showing the component due to the non-condensate (dashed). Param- 
eters used were T — 800nK, 9 million condensate atoms, and a 
density scale factor of 1.76 x 10~ 14 . 

Because the non-condensate density is taken to be given 
by a distribution with zero chemical potential, which means 
in part that the mean field interactions with the condensate 
are neglected (as now V^g(r) = Vt(t)), this allows the non- 
condensate distribution to be quite strongly peaked at small 
z, whereas the inclusion of this effect (as shown in Fig.[l7|) 
causes the distribution in this region to be suppressed. This 
could lead to an underestimate of the number of condensate 
atoms, which might be particularly significant at the initial 
stages of the condensate formation due to the low numbers 
present. Using the full MIT procedure with up to 5 extra pa- 
rameters must give a better fit, but this basic problem will still 
be present, since the qualitative behavior of the vapor in the 
condensate region will be very similar to this fit. 



21 



c 

CL 



c 
C 



o 
U 
— 

B, 




-400 



-200 



200 



400 





z (Hm) 

FIG. 20. Equilibrium data as for the previous figure, fitted by us- 
ing the semiclassical method proposed in this paper. Temperature 
and number of condensate atoms were taken to be the same as was 
used previously (T = 800nK and no = 9 x 10 6 ), however the 
density scale factor was changed to 4.72 x 10 -14 . 

The same data can be fitted using the distributions from 
equation (|84|). Figure |o| shows the same equilibrium data, 
fitted by using the same parameters but by using the semiclas- 
sical distribution — as before, we use the T and no determined 
by the MIT group using their procedure, so that the only free 
parameter is the scale factor. This spatial distribution can be 
seen to fit the experimental data much more closely. However 
this can still be improved using a different choice of T and 
no. The best fit was found with parameters of T = 900nK 
with 4 million condensate atoms and a density scale factor of 
5.4 x 10~ 14 , and is shown in Fig.|2l]. As a measure of how 
good the fit was, the sum of the square of the differences be- 
tween the data and the fitted curve was calculated. The best 
curve gave a value of 5.7 x 10 3 , the fit in Fig|2(] gave a value 
of 1.23 x 10 4 and the fit of Fig.pjlgave 1.76 x 10 4 . 
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FIG. 21. The best fit to the same data used in the previous two 
figures. The parameters used were T = 900nK, no = 4 x 10 6 and 
the density scale factor was 5.4 x 10~ 14 . 



The temperature determined by our method at 900K is not 
very different from the 800K determined by the MIT proce- 
dure, but the condensate number by our method is at 4 x 10 6 
less than half of the MIT value of 9 x 10 6 . Note however that 
the fit of Fig. ^ does not appear visually unacceptable, lead- 
ing to the conclusion that the best we can say is that the data 
do not constrain the condensate number very strongly. The 
reason for this probably lies in the choice of the least squares 
fitting procedure, which weights all data equally, and since 
there is not very much data in the condensate region, the con- 
densate number can be varied significantly without changing 
the goodness of fit greatly. 

In order to appreciate better the differences between the 
three fits, we have plotted all three fit functions on the same 
graph in Fig. 22 From this it is clear that the fits using the 
zero chemical potential Bose-Einstein Fig.[l2](a) and the more 
correct Fig.[22|(b) are much more similar to each other than to 
the best fit, Fig.^2|(c), in which the higher temperature makes 
for a much more pronounced peak in the vapor distribution, 
upon which a condensate peak is then superposed, giving a 
larger overall peak than would be expected from the conden- 
sate. From this figure, it is clear that the properties of the 
distribution in the transition region can be very deceptive. 

We should also note that the density scale factor does not 
change significantly from Fig.^ to Fig|2l]. Since this scale 
factor is a well defined physical quantity which should be de- 
termined more or less equally by all data, this is encouraging. 
By this we mean that, although the scale factor is of no par- 
ticular interest to the physics of Bose-Einstein condensation, 
it is nevertheless a well defined physical quantity whose value 
should be determinable. In contrast the fit of Fig. gives a 
very different scale factor, which is not unexpected, given that 
it does not pass through most of the data points. 
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FIG. 22. Comparison of the three fit functions: (a) Solid line, The 
zero chemical potential vapor fit as in Figjlj^; (b) Dotted line, Fit us- 
ing a nonzero chemi cal p otential vapor as in Fig.^; (c) Dashed line, 
The best fit as in Fig 
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FIG. 23. Comparison between experimental data from [|4|] and theoretical spatial density distributions, calculated using the semiclassical 
density distributions and the growth model described in this paper. The first frame shows the distribution before the cooling 'cut' below the 
critical temperature was performed. The parameters used w ere T — 800nK and no,/ = 9 x 10 6 . The (lower) red lines show the theoretical 
curves using initial condition (c) described in section [II A 6 In the second and third frames the (upper) green curves depict the results obtained 
using initial condition (d) which become essentially indistinguishable from the condition (c) lines after 40ms. 
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FIG. 24. Further theoretical spatial distributions as the condensate grows. These grap hs lead on from those in Fig.[23| and were obtained 
using the same parameters. The scale has been chosen so as to correspond to that in Fig]23. The last graph shows the distribution at equilibrium 
(400ms - solid line) along with the result at 160ms, the latest time for which the experimental data was available. Note the change of scale in 
the last graph. 
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D. Comparison of Spatial Distributions with Experiment 



With the aid of the spatial distributions calculated using 
equation (Q), comparisons with the MIT growth data are pos- 
sible. In figure 2 of [Q] the density profiles of the system are 
given for a single condensate growth. The profiles are one- 
dimensional slices through the centre of the system, in the z 
(long) axis. The fits to the equilibrium profile were compared 
in the previous section. From this data, MIT extracted val- 
ues of T = 800nK and n j = 9x 10 6 , and they found that 
the final equilibrium situation seemed to have been reached 
by 160ms [^2[]. The theoretical spatial distribution curves are 
compared to the experimental data in Fig.^3| The figure shows 
very good agreement with the experimental data during this 
160ms period. The frame at 40ms shows the least good agree- 
ment, however the agreement can be improved by the choice 
of different initial conditions as is shown in Fig.O. Using ini- 



tial condition (d), described in section |III A 6[ the agreement 
with experiment again becomes quite respectable. It is inter- 
esting to note that the spatial density distribution in this frame 
depends quite strongly on the initial conditions, which were 
found to have on ly a very small effect on the growth curve 
This is because the 40ms frame is taken 



III A 6 



in section 

at a time very close to the initiation of the fast growth of the 
condensate population, and the spatial distribution at this time 
is quite strongly dependent on exactly when the growth does 
start, since any difference in the initiation time creates a rela- 
tively large change in the occupation numbers given that they 
are still quite low at this time. 

The condensate growth curve which corres pon ds to the spa- 
tial distributions in Fig.^3] is shown in Fig 25 The figure 
shows that at 160ms (the time of the last graph in Fig.^3|) the 
growth of the condensate is theoretically only about halfway 
to completion. No experimental data for this particular growth 
event at times greater than 160ms was published in [ffl], al- 
though two further frames were taken at some later times and 
were found by the MIT group to contain the same number of 
atoms as for the 160ms profile [22|. 
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FIG. 25. The condensate growth curve corresponding to the situ- 
ation in Fig.^3[ The dashed line indicates the value at t = 160ms 
which corresponds to the final graph in Fig.^3| Note that the growth 
starts to become significant at about the t — 40ms frame. 



The spatial density distributions for times later than 160ms 
are shown in Fig.p_4|, for times up until 400ms at which time 
equilibrium has been reached in the theoretical model. The 
graphs in this figure show that the spatial distribution changes 
relatively little after 160ms compared to its earlier evolution. 
In particular the last graph in the figure compares the spatial 
distribution at equilibrium with that at 160ms, and it shows 
only a fairly small change. The peak height has definitely be- 
come greater, although not by a great deal more when one 
considers the size of the fluctuations evident in the experi- 
mental data of Fig.^3]. However, more importantly, the spatial 
extent of the condensate has not changed significantly after 
160ms. The condensate number was determined in the ma- 
jority of the MIT fi ts fr om the width of the condensate, as is 
given by equation (102). A well known difficulty [^9|] in us- 
ing this method is that, in the Thomas-Fermi approximation, 
the condensate occupation depends upon the fifth power of 
the condensate radius. Thus, for example, a 15% error in the 
measurement of the radius will result in the value for no be- 
ing incorrect by a factor of 2. When this is considered, along 
with the fact that the experimental data will show fluctuations, 
and given the comparison between the equilibrium and 160ms 
distributions, it seems perhaps possible (though in the opinion 
of the experimenters p^ ] very unlikely) that there was some 
evolution at times later than 160ms which was not able to be 
measured experimentally to sufficient accuracy. 

In conclusion therefore, more experimental data is really 
needed in order to make better comparisons, especially needed 
is data for times at which equilibrium has definitely been 
reached (at long t, although for very long times processes such 
as trap-loss may have to be considered). It would also be ben- 
eficial to fit the full two-dimensional data, rather than the one- 
dimensional slices shown previously. This would reduce the 
influence of fluctuations, and would lead to more weight being 
given to the measurements of the thermal cloud, which would 
increase the accuracy of any temperatures found. 



V. COMPARISON WITH OTHER THEORETICAL 
TREATMENTS 

Most of the other theoretical treatments have attempted to 
describe the formation of a condensate in a homogeneous, un- 
trapped situation. No real quantitative predictions have yet 
emerged from any of the work that has been performed on 
trapped dilute atomic gas BEC, and so comparisons unfortu- 
nately will have to be qualitative at best. 



A. Quantum Boltzmann Equation Approach 

One of the techniques used to describe condensate growth 
has been the quantum Boltzmann equation, which has been 
used by Snoke and Wolfe [j30|, Semikoz and Tkachev [|3l|], 
and Holland, Williams and Cooper [p0|], as well as forming 
the basis of the theory of the kinetic stages in the work of Ka- 
gan, Svistunov and Shlyapnikov [[32 33]. Although the theory 
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described in thispaper was developed from the Quantum Ki- 
netic theory [§-[ljJ], it turns out that essentially the same equa- 
tions may be obtained by modifying the quantum Boltzmann 
equation approach as follows: 

(i) The quantum Boltzmann equation in an ergodic form is 
used, a form similar to that used in 



df(e„) 8ma 2 w 2 



dt itTi 

/(e P )/(e,)(l + /(e m ))(l + /(e„)) 

-/(em)/(e„)(l + /(ep))(l + /(e,)) 

(103) 

where = gufi^k) is the number of particles with energy 

ejfe, e min = min(e m ,e„,e p ,e 9 ) and AE = e m + e n -e p -e q . 

(ii) The energy leve ls in the condensate band are modified as 
discussed in section III A 7 , in order to account for the mean 



field interactions due to the presence of the condensate. 

(iii) The levels in the non-condensate band are summed over 
and assumed to be time independent. This allows much larger, 
and realistic sized, systems to be modelled as opposed to the 
100-1000 atom systems typically simulated in previous at- 
tempts. 

(iv) Collisions between two particles which were both ini- 
tially in the condensate band were neglected, this is a valid 
approximation if the vast majority of particles are found in 
the non-condensate band. 

Using these modifications, and the rates for the scattering 
and gro wth processes found from quantum kinetic theory (see 
sec. II F ), the quantum Boltzmann equation will give rise to the 
set of differential equations ^% whose solutions provided the 
results in this paper. 

Of the above references, only the work of Holland, 
Williams and Cooper conducted any simulations for the 
growth of a trapped condensate. They found that their sim- 
ulations of condensate occupation number evolution behaved 
as no = «o,/(l — e~'/ T ), where r was a fitted parameter. A 
function of this form can be made to fit the results obtained 
by our model reasonably well, provided that an initiation time 
is allowed for, as was anticipated might be necessary in [|20|]. 
The same functional form was obtained using the quantum 
Boltzmann master equation approach in [|o|] by Jaksch et al. 

It should be noted that, although the quantum kinetic de- 
scription for the growth of the mean occupation numbers turns 
out to also be described by a modified quantum Boltzmann 
equation approach, the full quantum kinetic theory treats as- 
pects of condensate dynamics which are not accessible via the 
quantum Boltzmann equation. Such aspects include: the treat- 
ment of fluctuations, phase and phase decoherence, and the 
inclusion of Bogoliubov-like quasiparticle states. 



B. Comparison with Work of Kagan et al. 

Major theoretical work into the growth of condensate in 
recent years has been performed by Kagan, Svistunov, and 
Shlyapnikov. They divided the growth into three stages, the 
first of which was a kinetic stage described by Svistunov 
in [p2[]. Svistunov predicted the formation of a particle- 
flux wave in energy space during the initial stages, which 
transports particles towards lower energy states. The arrival 
of this wave at the lowest energy state at the critical time 
would give rise to an energy distribution function of the form 
f(E) oc E~ 7 / 6 . After the critical time this behaviour would 
be lost due to a particle-flux wave propagating to higher ener- 
gies. The simulations of Semikoz and Tkachev [gTJ showed 
this behaviour to some extent, although they found that the 
behaviour at the critical time was f(E) oc E" 1 - 24 . 

The work by Svistunov in [[32|] related to the case of homo- 
geneous systems, however he recently reworked his method- 
ology to consider a gas confined by a harmonic trapping po- 
tential. In this case he found [ p4[ | that the dependence at the 
critical time now tended towards f(E) oc E 1-5 / 3 . 
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FIG. 26. The growth of a condensate of 9 million atoms as pre- 
dicted by the model in this paper. Plotted is / (e m ), the population 
per individual level with energy Cmi and the energies of the levels 
as functions of time. The values of f(e m ) and e m are plotted on 
log scales. The lines almost parallel to the time axis are not lines 
of constant energy, but rather lines of constant level number, whose 
energies change with time. Not all levels are shown in order to make 
the behaviour legible. The solid black curve represents the energy of 
the condensate level log 10 [/ijv(no)] as a function of time. The bold 
line labelled C-C represents the critical time, at which the energy 
distribution has the form f(E) oc E~ 161 . The temperature of the 
system was 800nK. 

In Fig.^6] the results our model of condensate growth are 
shown in a somewhat different form. The occupation numbers 
for both the condensate level and excited states are plotted as 
a function of their energy and the time (note the logarithmic 
scales). From this graph several points can be noted, firstly 
the front corner shows quite how rapidly the initial conditions, 
however arbitrary, are smoothed out by scattering processes, 
and the discontinuous initial conditions rapidly approach a re- 
alistic distribution. The growth of the condensate is rather 
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small up to the point labelled as the critical time, after which 
the condensate grows enormously. The populations of the ex- 
cited states approach equilibrium very rapidly after the critical 
time, much more rapidly than the condensate level does. 

As the critical time is approached, the distribution ap- 
proaches a straight line as shown in figure |6| At the crit- 
ical time, when this distribution is linear with the logarithm 
of the energy, the energy dependence was found to be of the 
form f(E) oc E~ 1S1 , which is in good agreement with the 
prediction of Svistunov of E~ 5 / 3 . 

As far as estimates of timescales are concerned, Kagan, 
Svistunov and Shlyapnikov predicted that the evolution of the 
system up until initiation would occur on the timescale of the 
classical collision time to, that the initiation of the condensate 
growth would occur on a much faster timescale, and that the 
final growth would occur on the timescale necessary for the 
annihilation of vortices in homogeneous gases, and the decay 
of fluctuations in the phase of the condensate. 

In the treatment of Kagan, Svistunov and Shlyapnikov'the 
timescale for the first kinetic stage was postulated for a ho- 
mogeneous gas, however a comparison can still be made. A 
first estimate of tq = (anvx)' 1 can be obtained by using the 
classical value for vt, the mean thermal velocity in a gas, of 
Vt = \/2kT/m, and the cross section defined by a = 8wa 2 . 
The value for n, the mean density, will be taken as the density 
of the non-condensate particles in the centre of the trap given 
by 
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(see equation p8[) where G a {z) — J2^Li e qz /q a - Taking 
Yr(0) = 0, and using a temperature of 900nK, the collision 
time is tq = 27ms. This temperature corresponds to that of 
the growth in Fig.[| which shows that the time until initiation 
is of the order of 2tq. Thus our treatment does agree with the 
picture of Kagan et al. in that this stage occurs over the order 
of a few ro. 

A note about the collision times is needed here. In exam- 
ining the validity of the ergodic approximation Jaksch et al. 
found that it was valid only for quantities averaged over about 
10 collision times JlOj] . The above timescale tq is not the 
timescale over which collisions occur in the condensate sys- 
tem in reality. It is rather the classical collision time for a 
classical gas in equilibrium below the critical temperature but 
with no condensate present, and is obviously artificial. Once 
the condensate begins to form the density increases signifi- 
cantly, and the actual mean collision time was found in [ [l0| ] to 
be more than two orders of magnitude smaller than tq . Thus 
the ergodic assumption should still be valid for our treatment, 
even given the large value of tq. 

The timescales found by Kagan et al. for the second and 
third stages of evolution have only been determined for the 
case of a homogeneous gas, and so accurate comparisons with 
our model for these stages are not able to be performed. Our 
model does agree that the initiation stage occurs on a much 
faster timescale than the first kinetic stage. The presence of 
vortices has not been considered in our treatment, nor has 



any consideration been given to the phase fluctuations, and so 
comparisons cannot be made with the third stage of evolution 
in the description of Kagan et al. 



VI. CONCLUSIONS 

In this paper a model of the growth of a condensate, de- 
rived from quantum kinetic theory [[)- 1 1 1, is presented. The 
main improvements over the simpler description in Q| are the 
more accurate calculation of the growth rate factor W + (no), 
the consideration of the time dependence of the lower energy 
levels, and the inclusion of the scattering of particles between 
these levels. The modified W + factors have the greatest ef- 
fect, generally increasing the rate of growth by a factor of 
about 3, dependent on the exact parameters. The inclusion of 
the other levels and their scattering also leads to an increase 
in the rate of growth, mainly by reducing the amount of time 
taken for the initiation of the growth, that is the time before 
the stimulated growth processes, due to the Bose statistics of 
the system, become dominant. 

The model describes the evolution of time-dependent en- 
ergy levels in the lowest states, coupled with a time indepen- 
dent thermal bath of atoms occupying the higher energy lev- 
els. 

The results give growth curves whose shapes are approxi- 
mately those given by 



n 



for t < tj 

n QJ (l -e-C-**)/ 7 ") fori > U 



(105) 



where tj is some initiation time. This form agrees with the 
general form of the results of Holland et al. [ P-Oy and Jaksch 
et al. JlO| ] once an initiation time is allowed. The results also 
seem to be in qualitative agreement with features of the de- 
scription proposed by Kagan, Svistunov and Shylapnikov. 

The results are not very sensitive to the exact nature in 
which the mean-field effects on the lower levels were ac- 
counted for, so long as the energies of the very lowest lev- 
els were altered in a consistent fashion. The initial conditions 
used did not have a large effect on the growth curves, however 
they can be important when the spatial density profiles of the 
system are calculated for comparison with experiment. 

The evolution of the model depends upon approximations 
made for the rate factors T(T) and W+ + (no). The results ob- 
tained show that the growth is not very different as long as the 
actual value of T(T) falls within about a factor of 10 of the 
approximation, and as long as the actual values of W^ + {no) 
lies within a factor of 2 of W + (no). 

Overall the rates of growth predicted now agree very well 
with the growth rates measured from experimental data [Q]. 
However at lower temperatures the trend in growth rates 
shows some divergence, with the experimental rates becoming 
quicker whilst the predicted rates become slower. This may be 
a result of the substantial cooling necessary to achieve these 
temperatures, giving rise to a highly non-equilibrium system 
which is inadequately described by our model. 
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The extraction from experimental data of energy level pop- 
ulations is fraught with difficulties, leading to large uncertain- 
ties. It was found to be better to compare theoretical predic- 
tions of the evolution of the spatial density distribution of the 
system as the condensate to those experimentally measured, 
and the results in sect.0 appear to give good agreement with 
the spatial shapes given in the MIT data. They also show that, 
once a large condensate occupation has been reached, if the 
condensate occupation changes by a significant amount the 
corresponding change in the spatial density distribution of the 
system is relatively small. 

With the aid of more raw experimental data to fit, a deter- 
mination of exactly how accurate the predictions of this model 
are, conclusions as to the exact initial conditions present, and 
more accurate approximations of T(T) and W+ + (no) may be 
able to be reached. 

Further work which could be undertaken within the frame- 
work of this model includes: 

i) An accurate determination of the rate factors T(T) and 
Wj^ + (no) analytically, or at least finding constraints on their 
values by comparison with more experimental data. 

ii) An inclusion of some Bogoliubov phonon-like quasiparti- 
cle nature in the description of the lower energy levels, since 
all excited levels in this paper were treated as Hartree-Fock 
particle-like quaiparticles, which will be valid for most of the 
higher levels but not for the lower energy excitations. 

iii) A consideration of the fluctuations in occupation num- 
bers. These may be significant in determining the initiation 
time, which is when the occupation of the condensate level 
becomes large enough for the stimulated growth processes to 
take over. 

iv) As the model stands, the non-condensate band 'bath' of 
atoms is treated as being time independent. A major exten- 
sion of the model would be to include the dynamics of the 
non-condensate band in the evolution. Extension to include a 
time-dependent bath will be treated elsewhere [p5[]. 
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